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(54) Tide: MEASUREMENT OF GENE EXPRESSION PROFILES IN TOXICITY DETERMINATION 



(57) Abstract 



A method is provided for assessing the toxicity of a compound in a test organism by measuring gene expression profiles of selected 
tissues. Gene expression profiles arc measured by massively parallel signature sequencing of cDNA libraries constructed from mRNA 
extracted from the selected tissues. Gerte expression profiles provide extensive information on the effects of administering a compound to 
a test organism in both acute toxicity tests and in prolonged and chronic toxicity tests. 
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MEASUREMENT OF CfeNlE EXPRiESSION PROFILES 
' IN TOXICITY DETERMINATION 

'' ' ' * > . 

' Field of the Invention ' 



The invention relaties generally to inetliods' for detecting and monitoring 
phenotypic changes in in vitro and in vivo systems for assessing ahd/6f determining 
the toxicity of cheihicai coriripouiids, and riibre paiticidarly^the to a 

method for detecting and monitoring' changes in gene (expression patterns in in vitro 
1 0 and in vivo systems for detemiihihg the toxicity of drug candidates. ' ' 



BACKGROUND = > ^ ^ r - ■ 



The ability to rapidly and cdnveniehtl;j'' aisisess the tbxidty of hew compounds 

is extremely iniportant. Thbiisahds of hew compounds are synthesized every year, 

1 5 and many are introduced to the environment through the development of new 

commercial products and processes, often with littl^ knowledge of their ^ort term 
and long term health effects. In the development of new drugs, 'the cost of assessing 
the safely and efficacy of candidate conipounds is becoming astronomitai: It is 
^estimated that the phariiiaceuticai indiisti^ spends aii average of about 300 million 

20 dollars to bring a new'phanhaceutical compound to market, e.g: Biotech^ 

226-228 ( 1 995). A large fraction of these costs are due to the failure of candidate 
compounds in the later stages of the developmental process. That is, as the 
assessment of a candidate drug progresses from the identification of a compound as a 
drug candidate— for example, through relatively inexpensive binding assays or in vitro 

25 screening assays, to pharmacokinetic studies^ to toxicity studies, to efficacy studies in 
model systems, to preliminary clinical studies, and so on, the costs of the associated 
tests and analysies increases tremendously. Consequently, it may cost several tens of 
millions of dollars to determine that a once promising candidate comf)Ound possesses 
a side effect or cross reactivity that renders it commercially infeasible to develop 

30 further. A great challenge of pharmaceutical development is to remove from further 
consideration as early as possible those compounds that are likely to fail in the later 
stages of drug testing/ 

Drug development prograrrs are clearly structured with this objective in mind; 
however, rapidly escalating costs have created a need to develop even more stringent 

35 and less expensive screens in the early stages to identify false leads as soon as 

possible. Toxicity assessment is an area where such improvements may be made, for 
both drug development and for assessing the environmental, health, and safety effects 
of new compounds in general. 
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Typically the toxicity, of a coqipoimd is 4etenniqed by ,adjninistering the 
conipound to one or more species of t^t animal under controlled conditipns.and by 
monitoring the effects on a wide range of parameters. The parameters include such 
things as blood chemistry, ^^reight=gaifl or loss, a variety.of behavioral patterns, ^nuscle 
ic^evbody tonperatine, respiration rate, , lethalio^,and the like., which collectively - 
provide a measure of the state of health of the test animal. The degree of deviation of 
such parameters from their nonnal ranges, gives a measure of thjs toxicity of a. 
cpippountl. Such tests may be designed tp assess ^he acute,, prplonged..or chronic 
^Wf^. Wra^Pund, In general,.acute t^, inyolye admiiiistratipn of t^^ test 
chemical on one occasion. The period of oteervation of the test animals may be as 
short as a few hours, although.it is,usually at. least 24. hours and.in some cases -it may 
as loiig as a week.or.more.. In general. pro|pnged.tests involve administration of 
the jest chemical pri multiple occasions, Theiest chemical may be administered one 
or more times each day, irregularly as .when it is incorporated in the diet,,at sp«>cifiCi v.F 
times such as during pregnancy, or in some cases^regularly bpt only . atweekly." 
intervals. Also, in the prolonged test the experiment is usually conducted for not less 
than 90 days in the rat or mouse or a year in the dog. In contrast to the acute and 
prolonged types of test, the chronic toxicity tests are those iri which the test chemical 
- administered for a substantial portion of the iifetime of the test animal., In.the case 
20 of the mouse or, rat, this is a period of 2 to 3 years. In the case of .the, dog. it is for 5 to 
7 years. 

Significant costs are incurred in establishing and maiiitaining large cohorts of 
test animals for such assays, especially the larger animals in chronic toxicity assays. 
Moreover, because of species specific effects, passing such toxicity .tests does not 
25 ensure that a compound is free of toxic effects when u«:d in humans, Such tests do. 
however, provide a standardized set of information forjudging the safety of new . 
compounds, and they provide a database for giving preliminary assessments of related 
compounds. An important area for improving toxicity determination , would be the , 
identification of new observables which are predictive of the outcome of the 
30 expensive and tedious animal assays. 

In other medical fields, there has been significant interest in applying recent 
advances in biotechnology, particularly in DNA sequencing, to the identification and 

genes in healthy and diseased organisms, e.g. Adams 
et al. Science, 252: 1651-1656 (1991); Matsubaraet al. Gene, 135: 265-274 (1993); 

patent application, PCT/US95/01863. The objectives 
of such applications include increasing our knowledge of disease processes, 
identifying genes that play important roles in the disease process, and providing 
diagnostic and therapeutic approaches that exploit the expressed genes or their 
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products. Wliire such* approaches are attnictive/thbs^^ or even 

'sampled, sequencing of.expressed genes are still beset by the enormous effort : 
required: It is estimated that 30-3 5 ^thoiisanci diff^^ in a typical 

mammialiah tissue in ahy given' state, e.g; Ausubel et ai; Editors, Current Protocols, 
5 5:8.i-^:8.4 (John Wiley & Sd^^^ Determining the sequences of 

even a small sample of that number of gene products is -a major enterprise, requiring 
industrial-scale resources. 'Thus;' the^ro application of'massive sequencingfof 
expressed genes is stiirbe^^ . ? ; ; . ^ 

The iavailability of new iuisays for assess 
10 'as candidate drugs, thiat Voiild provide more compreherisrve* aind precisie information 
aboiit the state bf ifieialth of a test animal would be highly desirable. Such additional 
assays woiild'preferably be less expen^^^ 

ciirrerit teistirig jiroceduiies, and would at the 'same time provide enough infoimation to 
' make early judgmerits^^^^ 0 : > 

Summary of the Invention ' ' • j 
An object of the invention is to provide a new approach to toxicity assessment 

based on an examination of gene expression patterns, or profiles, in in vitro or in vivo 

''•testsystemis. '■■^-■■•■*. ^ . j*--' \' ■ i:- . 

20 Another object of the invention is to provide a database on which to base 

decisions concerning the toxicological properties of chemicals, particularly drug' 
candidates. ' ■ . < 

A further object of the invention is to jjrovide a method for analyzing gene 
expression patterns in selected tissues of test animals. 
25 A still fiirther object of the invention is to provide a system for identifying 

genes which are differentially expressed in response to exposure to a test compound. 
Another object of the invention is to provide a rapid and reliable method for 

correlating gene expression with short term and long term toxicity in test animals. 
if . ' . 

Another object of the invention is to identify genes whose expression is 

30 predictive of deleterious toxicity. "\ 

The invention achieves these and other objects by providing a method for 
massively pairallel signature sequencing of genes expressed in one or more selected 
tissues of an organism exposed to a test compound. An important feature of the 
invention is the application of novel DNA sorting and sequencing methodologies thai 

35 permit the formation of gene expression profiles for selected tissues by determining 
the sequence of portions of many thousands of different polynucleotides in parallel. 
Such profiles may be compared with those from tissues of control organisms at single 
or multiple time praints to identify expression patterns predictive of toxicity. 
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iKC of oligonucleotidie tags 
that are members of a minitnally CTpss-hybridlziog set pf oligonucleotides^ . The 
sequences pf oligonuclec^des the sequences of every other 

;,™5??^f ?f -^^set by atJeasttw^^ Thus, each member of such a set 

""^^^^ 4"I?Jex (or triplex) with the.complemem of any ptherm^raber with less 

„ Complement 

; , '^^Plements.; inasjcompri^ natu^ or.non,natural 

nucleotideanalogs. , Inferably, tag complemente.^^ solid.phase supports. 

, olig<?smicI^_de_tags when used with their coiresponding tag complements 
provide a. means of enhancing specificity of hybridization for sorting polynucleotides, 
such ascDNAs. ... 



. . ^ . pie polynucleotides to be sorted each haye.an oligonucleotide tag,attached, 
.,^uch tfi^t different^polynucleotides have different tags. .As explained more fully 
below, this,condition is acWeved by employing a repertoire of tags substantially 

5 greater than th^ population pf polynucleptides and by taking a sufficiently, small r^ . 
sample of tagged polyijucleotides from the full ensemble of tagged polynucleotides. 
After such sampling, when the populations of supports and polynucleotides are. mixed 
under conditions >yhich permit specific hybridization of the oligonucleotide tags with 
-their respective complements, identical polynucleotides sort onto particular beads or 

0 regions. The sorted populations of polynucleotides can then be.sequenced on the 
solid phase support by a "single-base" pr "b^e-by-base" sequencing methodology, as 
described more fully below. . 

In one aspect, the method pf the invention comprises the following steps: (a) 
administering the compound to a test organism; (b) exuacting a , population of mRNA 

5 molecules from each of one or more tissues of the test organism; (c) forming a 
separate population of cDNA molecules from each population of mRNA molecules 
extracted from the one or more tissues such that each cDNA molecule of the separate 
populations has an oligonucleotide tag attached, the oligonucleptide tags being 
selected from the same minimally cross-hybridizing set; (d) separately sampling each 

0 population of cDNA molecules such that substantially all different.cDNA molecules 
within a separate population have different oligonucleotide tags attached; (e) sorting 
the cDNA molecules of each separate population by specifically hybridizing the 
oligonucleotide tags with their respective complements, the respective complements 
being attached as uniform populations of substantially identical complements in 

> spatially discrete regions on one or more solid phase supports; (0 detennining the 
nucleotide sequence of a portion of each of the sorted cDNA molecules of each 
separate population to form a frequency distribution of expressed genes for each of 
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the one or more tissues; and (g) correlating the frequency distribution of expressed 
genes in each of the one or more tissues with the toxicity of the compound. 

An important aspect of the invention is the identification of genes whose 
^expression is iwredictive of the|oxicity of a conlpbund. Once ?ucfi genes are ' " ' 
identified/ they may be employed in conventional assays, such as TCverse transcriptase 
polymerase chain reaction (Rt-PCR) assays for gene expression"; 



Brief Description of the Drawings 
Figure 1 is a flow chart represen^^ an algorithin ifor generating 
10 minimally cross-hybridizing sets of oligonucleotides. hu 

Figure 2 diagrammatically illusirates an apparatus for carrying out 
polynucleotide sequencing in accordance^ with the invention. 



, ' . Definitions 
1 5 "Complement" or "tag complement" as used herein in reference to 

oligonucleotide tags refers to an oligonucleotide to which a oligonucleotide tag 
specifically hybridizes to form a perfectly matched duplex or triplex. In embodiments 
where specific hybridization results in a triplex, the oligonucleotide tag may be 
.selectefd to be either double stranded or single stranded. Thus, where triplexes are 
20 formed, the term "complement" is meant to encompass either a double stranded 

complement of a single stranded oligonucleotide tag or a single stranded complement 
of a double stranded oligonucleotide tag. 

r . . , T^^ .^5^ "oligonucleotide" as used herein includes linear oligomers of natural 
or modified monomers or linkages, including deoxyribonucleosides, ribonucleosides, 

25 anomeric forms thereof, peptide nucleic acids (PNAs), and the like, capable of 
specifically binding to a target polynucleotide by way of a regular pattern of 
monomer-to-monomer interactions, such as Watson-Crick type of base pairing, base 
stacking, Hoogsteen or reverse Hoogsteen types of base pairing, or the like. Usually 
monomers are linked by phosphodiester bonds or analogs thereof to form 

30 oligonucleotides ranging in size froni a few monomeric units, e.g. 3-4, to several tens 
of monomeric units. Whenever an oligonucleotide is represented by a sequence of 
letters, such as "ATGCCTG," it will be understood that the nucleotides are in 5'^3' 
order from left to right and that "A" denotes deoxyadenosine, "C" denotes 
deoxycytidine, "G" denotes deoxyguanosine, and "T" denotes thymidine, unless 

35 otherwise noted. Analogs of phosphodiester linkages include phosphorothioate, 
phosphorodithioate, phosphoranilidate, phosphoramidate, and the like. Usually 
oligonucleotides of the invention comprise the four natural nucleotides; however, they 
may also comprise non-natural nucleotide analogs. It is clear to those skilled in the 
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art when oligonucleotides living natural or non-natural nucleotides mav. be 
employed, e.g. where, proM l>y enzymes is ca|le4 for,. usually oligonucleotid^^^ 
consisting of natural nucleotides are required. . , 

oligonucleqtitle^sti^ up the duplex form a double §tranded, structure with 

one other such that every nucleotide in. each strand undergoes W^tsqn-Crick 
^^pP^j'^.58,r?i^ a nucleotide in the oth?r strapd. .The term.alsp comprehends the 
pairii)g of nucleoside analogs, sifch aa^.dt^x^sjne, nufjjj^idei?, wit|i.2rainlrippuri|» 
bases, and the like, .tlat^ma3/ te means j . 

that the triplex consists of a perfectly matched duplex and . a :third strand in which 
every nucleotide undergoes Hoogsteen or. reverse Hoogsteen association with.a . 
basepair of the perfectly matched .Conversely, a "mismatch" in a duplex 

between a tag ^d m oligoiiucleotide means that a pjur^.w in the 

duplex or triplex fails to undergo W^tson-Crick an^ and/or reverse , , ., 

15 Hoogsteen bonding! 

As used herein, "nucleoside" includes die natural nucleosides, including 2'- 
deoxy and 2'-hydro3cyl fprnis, e.g. as described in Kpmberg and Baker, DNA 
Replication, 2nd Ed. (Freeman, San Francisco,. 1 992). /^Analogs" in reference to 
. nucleosides includes synthetic nucleosides haying modified basej moieties and/or 
modified sugar moieties, e.g. described by.Scheit, Nucleotide Analogs (John Wiley, 
New York, 1980); Uhlman and Peyman,. Chemical Reyiews. 90: 543-584 (1990). or 
the like, with the only proviso that they are capablepf specific hybridization. Such 
analogs include synthetic nucleosides designed to enhance binding properties, reduce 
complexity, increase specificity, and the like. 

As used herein "sequence determination" or "determining a nucleotide 
sequence" in reference to polynucleotides includes deierniination of partial as well as 
fiiii sequence infoiroation of the pplynucleotide.. That is, the term includes sequence 
comparisons, fingerprinting, and like levels of information about a target 
polynucleotide, as well as the express identification and ordering of nucleosides, 
30 usually each nucleoside, in a target polynucleotide. The term also includes the 

determination of the identification, ordering, and locations of one, two, or three of the 
four types of nucleotides within a target polynucleotide. For example, in some 
embodiments sequence determination may be effected by identifying the ordering and 
locations of a single type of nucleotide, e.g. cytosines, within the target polynucleotide 
5 "CATCGC ..." so that its sequence is represented as a binary code, e.g. " 1 00 1 01 ..." for 
"C-(not C)-(not C)-C-(not C)-C ..." and the like. 
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As used herein, the term "complexity" in reference to a population of ^ 
polyiiucleoiides means the number of different species of molecule present in the 
population:' " • - --^^^^ V ^;^:v.:- ;v.^.. — 

As used herein, ttie terms "jgene expression profile," /and "gene expression \o 
5 paherh" which is iisecl equivalehily, means a frequency distribution of sequences of ' i 
portions of cDNA ixiolecules sampled froni a population of tag-cDN A conjugates. 
Generally, the portions of sequence are sufficiently long to uniquely identify the 
cDNA from which the portion arose. Preferably, the total number of sequences rM^;. 
detennihed is iai least ' 1 OOt); more preferably, the total niimber of sequences ^ 
10 determined in a gene expression profi . r 

As used herein, "test organism" means any in vitro or in viv6 system which ■ 
provides measureable responses to exposiute to test compK)uhds. Typically^ test 
orjganisms may be mammalian cell cultures, particularly of specific tissues, siichias 
hepatocyteis, neurons, kidney -cells, colony forming cells, or the like, or test organisms 
15 may be whole animals, such as rats, mice,' hianlsters, guinea pigs; dogs, cats, rabbits, 
pigs, monkeys, and the like. 

Detailed Description of the Invention 
The invention provides a method for determining the toxicity of a compound 

20 by analyzing changes in the gene expression profiles in selected tissues of test 
organisms exposed to the compound. The invention also provides a method of 
identifying toxicity markers consisting of individual genes or a group of genes that is 
expressed acutely and which is correlated with prolonged or chronic toxicity, or 
suggests that the compound will have an undesirable cross reactivit>'. Gene 

25 expression profiles are generated by sequencing portions of cDN A molecules 
construction from mRNA extracted from tissues of test organisms exposed to the 
compound being tested. As used herein, the term "tissue" is employed with its usual 
medical or biological meaning, except that in reference to an in vitro test system, such 
as a ceil culture, it simply means a sample from the culture. Gene expression profiles 

30 derived from test organisms are compared to gene expression profiles derived from 
control organisms to determine the genes which are differentially expressed in the test 
organism because of exposure to the compound being tested. In both cases, the 
sequence information of the gene expression profiles is obtained by massively parallel 
signature sequencing of cDNAs, which is implemented in steps (c) through (0 of the 

35 above method. 

Toxicity Assessment 
Procedures for designing and conducting toxicity tests in in vitro and in vivo 
systems is well known, and is described in many texts on the subject, such as Loomis 
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et al. Loomis's Esstentials of Toxicology, 4th Ed. (Academic Press, New York, 1996); 
Echobichon, The Basics of Toxicity Testing (CRC Press, Boca Raton. 1992); Frazier. 
editor, In Vit^q JoMcity . Testing (MareelDekker, New.York.:.1992); and the like. . 
Jn toxicity, te^ting,^twp,grpups of test one 
^ 8^9"P; ?ery^ as a .control and the other group, receives the, test compound: in a single 
dose (for acute toxicity tests) or a regimen of doses (for prolonged or chrdnic toxicity 
;?'P<^?>W.?!»ost.<»ses, %extract^^ of tissue as called for.in. the method of the 
invemion require? sacrificing the test animal,. both the control group and the group . . 
receiving conipound;mustbe large enough to, permit removal ofanimals, for sampling 
1 0 tissues,, if it is desired to observe the dynamics of gene expression through'the 
;, duration.of m experiment. ^. , ^, r^, :•.•,.■,,,},,,,:,,. .ii.-,,, : 

,. - in setting «p;a toxicity study^extensive guidance.is provided in the literature 
for selecting .the appropriate test organism for the compound being tested, route of 
administration, dose ranges, and the lijke. Water or. physiological-saline (0.9% NaCl 
1 5 in water) is the solute of choice for. the, te^t, compound since these solvents permit 
. administration by a variety of routes, . When this is not,possible because of solubility 
limitations, it is necessaiy to resort to the use of vegetable oils such as com oil or ; 
even organic solvents, of which propylene glycol is commonly used. Whenever 
-possible the use of suspension of emulsion .should be avoided.except for oral 
20 adminisu^tion. Regardless of the route of administration, the volume;required to 

administer a gi ven dose is limited by the size of the animal that is used. It is desirable 
to keep the volume of each dose uniform within and between groups of animals. . , 
When rates or mice are used the. volume administered by.the oral route should not 
exceed 0.005 ml per gram of animal., Eyen when aqueous or physiological saline 
25 solutions are used for parenteral injection the volumes that are tolerated are limited. 

although such solutions are ordinarily thought of as being innocuous. The 
intravenous LD50 of distilled water in the mouse is approximately. 0.044 ml per gram 

and that of isotonic saline is 0.068 ml per gram of mouse. 

When a compound is to be administered by inhalation, special techniques for 
30 generating test atmospheres are necessaiy. Dose estimation becomes very 

complicated. The methods usually involve aerosolization or nebulization of fluids 
containing the compound. If the agent to be tested is a. fluid that has an appreciable 
vapor pressure, it may be administered by passing air Uirough the solution under 
controlled temperature conditions. Under these conditions, dose is estimated from tiie 
volume of air inhaled per unit time, the temperature of tiie solution, and the vapor 
pressure ofUie agent involved. Gases are metered from reservoirs. When particles of 
a solution are to be administered, unless the particle size is less tiian about 2 ^m the 
particles will not reach the terminal alveolar sacs in the lungs. A variety of 
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apparatuses and chambers are available to perforin studies for detecting effects of ' 
irritant or othier toxic endpoints when they are administered by inhalation/ The * 
preferred method of administering an agent to animals is via the oral route, either by*: 
intubation or by incorpoi-ating the agent in the feed. ^ ' " ' ^ - - - - i y 

5 . Preferably, in (iesigriing a toxicity assessment, two or more spiscies should be 
employed that handle the test compound as similarly to man as possible in terrhs of 
metabolism, absorption, excretion, tissue storage, and thelike'. Preferably, multiple • ' 
doses or regimens at different concentrations should be employed to establish a dose* 
response relationship with respect to toxic effects. And preferably, the route of 
1 0 administration to the test animal should be the same as, or as siniilar zs possible to; 
the route of administration of the compound to rnan^ Effects obt^iied by one route of 
administraition to test animals are hot a priori applicable to effects by another route of 
administration to man. For example, food additives foi- man should be tested by ' ' 
admixture of the rhaterial in the diet 6f the test animals. ' " * * ' ^ ^ 

' ' • . ,^ • : I * I ; 1 • , - 1 . I 

15 Acute toxicity tests consist of adniinistering a compound to test organisms on ' 

one occasion. The purjwse of such test is to deterniine the symptomotology^^ - 
consequent to administration of the compound and to determine the degree of lethality 
of the compound. The initial procedure is to perform a series of rarige-fmdihg doses 
-of the compound in a single species. This necessitates selection of a route of ' ^ 

20 administration, preparation of the compound in a form suitable for iadministration by 
the selected route, and selection of an appropriate species. Preferably, initial acute 
toxicity studies are performed on either rats or mice because of their low cost, their 
availability, and the availability of abundant toxicologic reference data on these 
species. Prolonged toxicity tests consist of administering a compound to test 

25 organisms repeatedly, usually on a daily basis, over a period of 3 to 4 months. Two 
practical factors are encountered that place constraints on the design of such tests: 
First, the available routes of administration are limited because the route selected 
must be suitable for repeated administration without inducing harmful effects. And 
second, blood,, urine, and perhaps other samples, should be taken repeatedly without 

30 inducing significant harm to the test animals. Preferably, in the method of the 
invention the gene expression profiles are obtained in conjunction with the 
measurement 6f the traditional toxicologic parameters, such as listed in the table . 
below: • - 



-9 



wo 97/13877 



PCT/US96/16a42 



Hematoloj 
erythrocyte count 
total leukocyte count 
differential leukocyte 
count 

hematocrit ' 
hemoglobin 



Blood C hemist 
sodium 
potassium i ; 
chloride 



'J' 



calcium 
carbon dioxide 
serum glutamine-pyruvate 
transaminase 

serum glutamin-oxalacetic 
transaminase' • 
serum protein 
electrophoresisi y - 
blood sugar 
blood urea nitrogen 
total serum protein 
serum albumin 
total serum bilirubin 



Urine Ana lyses 

pH 

specific gravity 
total protein 

sediment 
glucose, 
ketones 

bilirubin 



t ! 1 



Oligonucleo tide Tags and Tag'Cbmp tements: ' > 
Oligonucleotide tags are members of a miiiimally cross-hybridizing set of 
oligonucleotides. The sequences of oligonucleotides of such a set differ from the 
sequences of every other member of the same set by at least two nucleotides. Thus, 
each member of such a set cannot form i duplex (or triplex) with the compiemem of 
any other member with less than two mismatches. Complements of oligonucleotide 
tags, referred to herein as "tag complements," may comprise natural nucleotides or 
non-natural nucleotide analogs. Preferably, tag complement^ are attached! to solid 
phase supports. Such oligonucleotide tags when used with their corresponding tag 
complements provide a means of enhancing specificity of hybridization for sorting, 
tracking, or labeling molecules, especially polynucleotides. 

Minimally cross-hybridizing sets of oligonucleotide tags and tag complements 
may be synthesized either combinatorially or individually depending on the size of the 
set desired and the degree to which cross-hybridization is sought to be minimized (or 
stated another way. the degree to which specificity is sought to be enhanced). For 
example, a minimally cross-hybridizing set may consist of a set of individually 
synthesized 10-mer sequences that differ from each other by at least 4 nucleotides, 
such set having a maximum size of 332 (when composed of 3 kinds of nucleotides 
and counted using a computer program such as disclosed in Appendix Ic). 
Alternatively, a minimally cross-hybridizing set of oligonucleoUde tags may also be 
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assembled combinatorially from subunits which themselves are selected from^^^ 
minimally cross-hybridizing set, iFor example, a set of minimally cross-hybridizing 
12-mers differing/from one another by at least three nucleotides may be synthesized 
by assembling 3 subimits selected from a set of minimally cross-hybridizing 4^mers 
5 that-each differ froni one another by three nucleotides. Such aii embodiment gives a 
maximally sized setof 9^, or 729, 4 2-mers. The number 9 is number of ' 
olijgonucleotides listed by the computer prograxh iof Appendixla, which assunies, as 
with the 10-mers, that only 3 of the 4 different types of nucleotides are used. The set 
is described as ''maximal'' because the computer programs of Appendices la-c provide 

1 0 the largest set for a given input (e.g. length, composition, difference m hiimber of 
.nucleotides between members). Additional minimally cross-hybridizing sets may be 
formed from subsets of such calculated sets. 

Oligonucleotide tags may be single stranded and be desigrieo for specific 
hybridization to single stranded tag complements by duplex formation or for specific 

1 5 hybridization to double stranded tag complements by triplex formation! 

Oligonucleotide tags may also be double stranded and be designed for specific 
hybridization to single stranded tag complements by tnplex fbrmatiori. 

When synthesized combinatorially, an oligonucleotide tag preferably consists 
.of a plurality of subunits, each subunit consisting of an oligonucleotide of ^to 9 

20 nucleotides in length wherein each subunit is selected from the same minimally cross- 
hybridizing set. In such embodiments, the number of oligonucleoiide tags available 
depends on the number of subunits per tag and on the length of the suhiiinits. The 
number is generally much less than the number of all possible sequences the length of 
the tag, which for a tag n nucleotides long would be 4"^. : ^ - 

25 Complements of oligonucleotide tags attached to a solid phase support are 

i V - A » * - t ^ , . , ■ I ^ 

,used to sort polynucleotides from a mixture of polynucleotides each containing a tag. 
Complements of the oligonucleotide tags are synthesized on the surface of a solid 
phase support, such as a microscopic bead or a specific location on an array of 
synthesis locations on a single support, such that populations of identical sequences 

30 are produced in specific regions. That is, the surface of each support, in the case of a 
bead, or of each region, in the case of an array, is derivatized by only one type of 
complement which has a particular sequence. The population of such beads or regions 
contains a repertoire of complements with distinct sequences. As used herein in 
reference to oligonucleotide tags and tag complements, the term ^'repertoire" means 

35 the set of minimally cross-hybridizing set of oligonucleotides that make up the tags in 
a particular embodiment or the corresponding set of tag complements. 

The polynucleotides to be sorted each have an oligonucleotide tag attached, 
such that different polynucleotides have different tags. As explained more fully 
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15 



below, this cpndition is achieved by employing a repertoife of tags substantially 
*® population of ppjiynucleotide^.and A)y.taking a! sufficiently small 
sample of tagged pplynudeotides.frpm the full enseinble of laggedipolynucleotides 
After, such sampling, >yheri the populations qf supports and polynucleotides are mixed 
under concUtions which i^nmts^^ 

their respective complements, identic^ sort Onto.particular beads or 

regions. • • . 

. ; .Theiiucleotide sequences of oligonuc|«)tides of a minimally: cross-hybridizing 
set are conveniently eniimerated by ,simple,.computer progiams; such as those . 
exemplified .^)y prografps whose source codes are listed in Appendices la and lb 

cpinputes.all minimally cross-hybridizing sets having 
4-mer subimits composed pf three kinds of nucleotides, Program, tagN^of Appendix 
lb enumerates longer oligonucleotides of a ^^^^ 

algorithms and,computer programs are readily written for listing oligonucleotides of. 
minimally cross-hybridizing sets for, any embodiment of the invention. Table r below 
provides guidance as to the size of sets of minimally cross-hybridizing:);, n , .u,', 
oligonucleotides for th^ indicated lengths and number of nucleotide differences. The 
above computer programs were used to generate the numbers. . ; • 



1 . .11 



20 



1 k ^ * 



Table I 



Oligonucleotid 
c 

. ..Word 
Lcneth 



Nucleotide 
; Difference 
between 
Oligonucleotides 
of Minimally 
Cross- 
Hybridizing Set 



I > 1 



Maximal Size 
ofMinimally 

Cross- 
Hybridizing 

1 1 'Set 



'Size of 
Rcpenoire 
with Four 

Words 



. Size of 
Repenoire with 
Five Words 



4 


■ -3 . 


' ' 9 


6561 


5.90 X lO'* 


6 


. 3 . 


27 ' 


5.3x10*'' 


1.43 X 10^ 


7 


4 


27 


5 J X lo' 


1.43 X 10^ 


7 ' 


5 


8^ 


4096 


3.28 X I O"* 


8 




m 


1.30 X 10^ 


2.48 X lo" 


8 


4 


62 


1.48 X 10^ 


9.16 X 10* 


8 


5 


18 


1.05 X 10^ 


1.89 X 10* 


9 


5 - . 


- 39 


2-31 X 10^ 


9.02 X 10^ 


10 


5 


332 


1.21 X lo'** 




10 


6 


28 


6.15 X 10^ 


1.72 X 10^ 


II 


5 


J87 






18 


6 


«25000 
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For some embodiments of the invention, where extremely large repertoires of 
tags are not required, dligonuclebtide tags of a mihinially cross-hybriidizirig set may 
be separately synthesized, $ets, containing seyeral.hundred to several .thousands, or 
5 even several tens of thousands, of oligonucleotides may be synthesized directly by a 
variety.pf .parallel;synthesis.apprQaches^ e.;g. as disclosed in Frank ei al, U.S. patent 
4,689,405;'Frank et al. Nucleic Acids Research, 11: 4365-4377 (1983); Matson et al. 
Anal. Biochem., 224: JlO-1 16 (1995); Fodoret al. International application 
P!GT/US93/041,45; Pease etal, Proc. Natl.,Acad..Sci.,;91 : 5022-5026 (1994); 
1 0 , Southern et ,al, J. Biotechnolpgj;, 3,5 : 2 1.7-22.7 (.1 994), Brennan, International 

. application PqT/US94/05896; Lashkari et al, Proc. Natl. Acad.' Sci., 92: 7912-7915 
(1995); or the like. y j,,, , . 

Preferably, oligonucleotide tags of the invention are sjiithesized . 
combinatorially out of subunits between three and six nucleotides in length and 

1 5 .selected from the same. minimally cross-hybridizing set. For oligonucletides in this 
range, the: members of such sets niay be enumerated by computer programs based on 
. the algorithm of .Fig. 1 . , : . . , 

. The algorithm of Fig Lis implemented by first defming, the characteristics of 
the subunits of the minimally cross-hybridizing set, i.e. length, number of base^ 

20 differences between members, and composition, e.g. do they consist of two, three, or 
.four kinds of bases. A table M^, n=l, is generated (100) that consists of all possible 
sequences of a given length and composition. An initial subunit S i is selected and 
. compared ( 1 20) with successive subunits S\ for i=nt 1 to the end of the table. 
Whenever a successive subunit has the required number of mismatches to be a 

25 member of the minimally cross-hybridizing set, it is saved in a nevy table Mn+ ] (125), 
that also contains subunits previously selected in prior passes through step 120. For 
example, in the first set of comparisons, M2 will contain Si ; in the second set of 
comparisons, M3 will contain Sj and S2; in the third set of comparisons, M4 will 
contain Si , S2, and S3; and so on. Similarly, comparisons in table Mj will be 

30 between Sj and all successive subunits in Mj. Note that each successive table M^+i 
is smaller than its predecessors as subunits are eliminated in successive passes 
through step 130. After every subunit of table has been compared (140) the old 
table is replaced by the new table Mn+i , and the next round of comparisons are 
begun. The process stops ( 1 60) when a table is reached that contains no 

35 successive subunits to compare to the selected subunit Sj, i.e. Mn=Mn+i . 

Preferably, minimally cross-hybridizing sets comprise subunits that make 
approximately equivalent contributions to duplex stability as every other subunit in 
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10 



15 



20 



25 



30 



35 



the s?t. In this, way, the stability of perfi?ctly rpa^ched duplexes between eveiy subunit 
and its comjjlement is approximateiy equal., Guidance for selecting such sets is 
Ki?^.!^;!!!!*!-!^ for selecting optimal PGR primers, and calculatmg 

dMplex st^iiiues.^ et al. Nucleic Acids Research, 17:,8$43-8551 (1989>. 

and 18: 6409-6412 (199p); &^^^ Natl. Acad, Sci., 83: 3746-3750, 

0986); Weitmur, Crit.Rey. Biochem. Mol. Biol., 26: 227-259 (1991);and the like. 
For shorter tags, e.g. about 3^ I?ss, .Ae^dgoritlunidescribed by Rychlik 

f H P'^^ff ®**V^i^°'"''W:|?SS' e.g. about 3.0.35 riucleotides or greater, 
an algorithm d^scjo^^^ Suggs et al, pages 683-6M in Brown, editor, I.CN-UCLA 
Syrnp. Dev. Biol., Vp^^^^ 

employed. Clearly, the are niany approaches ayailajble to one skilled in the art for 
designing s^s,ofmijumally cross:hybridi;ang subunits within the scope of the 

For example, to minimize the' affects of different base.stacking energies of 
tenminal nucleotides when subuiyts are .assembled, subunit? -may beiproyided.^^^ . 
have the same teminal nucleotides. In tliis- way,, when sul)units,are linked, the sum of 
the base-stacking energies of ail the adjoining terminal nucleotides will be the same, 
thereby reducing or ehininatii?g vari^ ,. 
_ A "word" of terminal nucleotides, shown in italic below, may. also, be added to 
each end of a tag.so that a perfect m^tch is always formed between it and a similar 
terminal "word" on any other tog complement. Such, an augmented tag woujd haye 
the form: 



w 


w, 


W2 . ..... Wk., . 




w . 


w . 




Wj' ... w.v 




w 



where the primed W's indicate complements. With ends of tags always forming 
perfectly matched duplexes,, all mismatched words will be internal mismatches 
thereby reducing the stability of tag-complement duplexes that otherwise would have 
mismatched words at their ends. It is well known that duplexes with internal 
mismatches are significantly less stable than duplexes with the same mismatch at a 
terminus. 

A preferred embodiment of minimally cross-hybridizing sets are those whose 
subunits are made up of three of the four natural nucleotides. As will be discussed 
more fully below, the absence of one type of nucleotide in the oligonucleotide tags 
permits target polynucleotides to be loaded onto solid phase supports by use of the 
5'-^3' exonuclease activity of a DNA polymerase. The following is an exemplary 
minimally cross-hybridizing set of subunits each comprising four nucleotides selected 
from the group consisting of A, G, and T: 
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5 Table II 

Word; W| W2 W3 W4 

Sequence: GATT TGAT TAGA TTTG 

Sequence: GTAA AGTA ATGT AAAG 

I'f '■ V, • ■ 

i . 

In this set, each member would form a duplex having three mismatched bases with. 
10 the complement of every. other member. - > * . ■ . . 

Further exemplary minimally cross-hybridizing sets are listed below in Table 
III Clearly, additional sets can be generated by substituting different groups of 
'nucleotides, or by using subsets of known minimally cross-hybridizing sets. 

15 . 

Table III 

Exemplary Minimally Cross-Hvbridizing Sets of 4-mer Subunits 



Set 1 


Set 2 


Set 3 


Set 4 


Set 5 


Set 6 


CATT 


ACCC 


AAAC 


AAhG 


AACA 


AACG 


CTAA 


AGGG 


ACCA 


ACCA 


ACAC 


ACAA 


TCAT . 


CACG 


AGGG 


AGGC 


AGGG 


AGGC 


ACTA 


CCGA 


CACG 


CACC 


' CAAG 


CAAC- 


TACA 


CGAC 


CCGC 


CCGG 


CCGC 


CCGG 


TTTC 


.GAGC 


CGAA 


CGAA 


CGCA 


CGCA 


ATCT 


GCAG 


GAGA 


GAGA 


GAGA 


GAGA 


AAAC 


GGCA 


GCAG 


GCAC 


GCCG . 


• ■ GCCC 






GGCC 


GGCG 


GGAC 


GGAG 
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Set / Set 8 Set 9 Set 10 Set 11 Set 12 

AAGA,-..,- •AAGC.-f,., .: ,,.AAGG,. , ACAQ.vm- c:|AC,CG...^v.....~ACGA~ 

^.9AC...^_. , ACAA • . ACAA , AACA , . . AAAA. AAAC 

AGCG - AGCG- ■ AGCC •• AG^t" '-" ^AGGc' - "'- AGCG'"- 

GAAG.m - ■ CAAGU - CAAG. - •> . CAAG .;:; i-..-;GACC 'J:v,:-,. T-CAGA.:r, -v. 

9PPA..-./ ......:CC.CC,, ,. . CCC.G, ,GCGA,. . .. ... ...GCGA • CCAG 

CGGC _ CGGA ■ CGGA ' ■CGCG ' ■ CGAG -CGGC ' 



iSAGA'- -.^GAGA - 'GAGG'?-'' ' • ''GAGGi"* • '-'i- GAGG'"- ■ •' 
GGGG, . ; . ,GGGG, ,., -.GCGe. ., .iGCCC,,.if,,.,..GCAG , GCCC,,. ,. 



d^^, .... r°?^F.-.r ^^AG , . GGAA GGCA GGAA 



oligonucleotide tiaig^^^^^ theif complements are'- ' ' 

conveniently srynthesize^ on ail autoiriatibd DNX synthesizfeK ci'iV an Applied ''''' 
Biosy^tems,' toe . (Foster Ci^; CaKfomiaV model 592' or ^94 DNA/RNA ^Syrithesizer: ' 
usih^ 'standard chemistries, sucfi aS phosphbi^Midite chemistry^'e;^/ disclosed in A 
following reiferehces-'iBeaucage and Iyer, Tetrahedr6h,'48: 2223-231 l"(i 992);' Molko 
et ii:'U.S: patent 4,980,460; Koster et al,' U.S. patent 4;725;677- Cariithers et al, U.S. 
patents 4,415,732; 4;458,066; arid 4,973,679; and the like. Altemativfe chemistries 
e.g. resiiltihg in non-natural backbone groups;' such as phosphorbthioke; ' ? ' 
phbsphbramidate, and the like, niay also be employed provided that the resiihing 
.oligonucleotides are capable of specifib hybridization. In some enibodimems, tags ' 
may comprise naturally occurring nucleotides' that permit jjrocessihg or manipulatibn 
by enzymes, while the eorrest)bnding tag complements rnay comprise hon-natiiral * 
nucleotide analogs, such as peptide nucleic acids, or like compounds, that promote the 
formation of more stable duplexes during sorting. ' 

When microparticles are used as supports, repertoires 6f oligonucleotide tags 
and tag cornplements may be generated by subunit-wise syntliesis via"split and mix" 
techniques, e.g. as disclosed iri ShorUe et al. International patent application ^ 
PCT/US93/0341 8 or Lyttle et al, Biotechniques, 19: 274-280 (1995). Briefly, the 
basic unit of the synthesis is a subunit of the oligonucleotide tag. Preferably, 
phosphoramidite chemistry is used and 3' phosphoramidite oligonucleotides are 
prepared for each subunit in a minimally cross-hybridizing set, e.g. for the set first 
listed above, there would be eight '4-mer 3'-phosphoramidites. Synthesis proceeds as 
disclosed by Shortle et al or in direct analogy with the techniques employed to 
generate diverse oligonucleotide libraries using nucleosidic monomers, e.g. as 
disclosed in Teleniiis et al. Genomics, 13: 718-725 (1992); Welsh et al, Nucleic Acids 
Research, 19: 5275-5279 (1991); Grothues et al, Nucleic Acids Research, 21 : 1321- 
1322 (1993); Hartley, European patent application 90304496.4; Lam et al, Nature. 
354: 82-84 (1991); Zuckerman et al, Int. J. Pept. Protein Research, 40: 498-507 
(1992); and the like. Generally, these techniques simply call for the application of 
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mixtures of tlie activated monomers to the growing oligonucleotide during the 
coupling steps. Preferably, oligonuclebtide tags and tag complements are synthesized 
on a DNA syrithesizer'havirig a num^^^^^ which is greater thanor 

equal to the number of different kinds of words iised in th€l construction of the tags. 
5 That iC preferably there is a synthesis chamber corresponding to each type of word. 
In this embodiment, words are aidded nucleotide-by-nucleotide, such 'that if a word 
consists of five nucleotides there i^eTive^mbrionieir cou^ each synthesis ' 

chamber. After a word is conipletely synthesized,' the s>^thesis^ sup removed 
from the c^lambersV mixed; arid redistiibuted batkto'the chambeii for^^ 
10 of word addition.' This latter embodiriient vSkki adVaritage of the high cbupling yields 



(5 



of monomer additiohre.g. in pliospKbraihi^^ ^ . u : \ 



Double stranded forms of tags may be made by separately synthesizing the 
complementary sWahds followed by' mixing lind^^ permit duplex 

' formation. Alternatively, double stranded tajgs may be formed by. first synAesizing a 
15 single stranded repertoire linked to a knb\wi 6Iigohu^^^ thatserves as a 

primer binding site.* The second strand is then syiithesized by combining the single 
stranded repertoire with a primer and extending with a polymerase: This latter 
approach is described in Oliphant et al, Gene, 44: 1 77-1 83 (1 986). Such duplex lags 
. may then be inserted into cloning vectors along with target polynucleotides for sorting 
20 and manipulation of the target polynucleotide in accordance with the invention; 

When tag complements are employed that are made up of nucleotides that 
have enhanced binding characteristics, such zs PNAs or oligonucleotide N3'->P5' 
phosphoramidates, sorting can be implemented through the fonriation of D-loops 
between tags comprising natural nucleotides and their PN A or phosphoramidate 
25 complements, as an alternative to the "stripping" reaction ennploying the 3'r->5' 
exonuclease activity of a DNA polymerase to render a tag single stranded. 

Oligonucleotide tags of the invention may range in length froni 12 to 60 
nucleotides or basepairs. Preferably, oligonucleotide tags range in length from 1 8 to 
40 nucleotides or basepairs. More preferably, oligonucleotide tags range in length 
50 from 25 to 40 nucleotides or basepairs. In terms of preferred and more preferred 
numbers of subunits, these ranges may be expressed as follows: 

Table IV 

Nimibers of Subunits in Tags in Preferred Embodiments 



Monomers 

in Subunit Nucleotides in Oligonucleotide Tag 

(12-60) (18-40) (25-40) 
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3 4-20subimits 6-13 subunits . 8-13 subunits 

! ; J:4i: V L.fcin . V3..15 subunits' ' ■ 4^10 subunits ' '6-10 siiijumts ' ' '' " 

2^:10 Siib\irtits ^ 3^muriite ^' ' 4^^ 

Most preferably, oligonucleotide tags are single stranded and specific hybridization 
occurs via Watson-Crick pairing with a tag complement. 

Preferably, repertoires of single stranded oligonucleotide tags of die invention 
contain at least 100 members; more preferably, repertoires of such tags contain at 
least 1000 members; and m6« prefCTablyrrepdioirW df such t^^^ contain at least 
10,000 members. r-yi.S i- ., "i-^-'-^/kilr^.'-- 

•••• -i'^ Triblex Taps '' ^-■■^l-K-^-l^ '■^^^'--^■■^ ' ' - 

In embodiments where" specific hybridizaiion-ocicurs via triplex fomiation, 
coding of tag sequences follows the same principles as for duplex^foraiing tags; 
however, there are further constraints on the selection of subiihit'sequencbs. 
Generally, third strand associaiioh via Hoogsteen type' of binding is most stoble along 
homopyrimidine-hbmopurine tracks in a double shanded target. ' Usu^ly, base triplets 
form in T-'A*T or C-G*G motifs (where indicates Watsbn-Crick paiiihg and 
indicates Hoogsteen type of binding); however, 'other 'niotifs afe also possible. For 
example, Hoogsteen base pairing permits parallel and antiparallel orientations 
between the third strand (the Hoogsteen strand) aiid the purine-rich strand of the ' 
duplex to which the third strand binds, depending on conditions and the composition 
of the strands. There is extensive guidance in the literature for selecting appropriate 
sequences, orientation, cbnditioris, nucleoside type (e.g. whether ribose or 
deoxyribose nucleosides are employed), base modifications (e.g. methylated cytosine, 
and the like) in order to maximize, or otherwise regulate, triplex stability as desired in 
particular embodiments, e.g. Roberts et al, Proc. Natl. Acad. Sci., 88: 9397-9401 
(1991); Roberts et al, Science, 258: 1463-1466 (1992); Roberts et al, Proc. Natl. 
Acad. Sci.. 93: 4320-4325 (1 996); Distefano et al, Proc. NaU. Acad. Sci., 90: 1 1 79- 
1 183 (1993); Mergny et al. Biochemistry, 30: 9791-9798 (1991); Cheng et al, J. Am. 
Chem. Soc, 1 14: 4465-4474 (1992); Beal and Dervan, Nucleic Acids Research, 20: 
2773-2776 (1992); Beal and Dervan, J. Am. Chem. Soc, 1 14: 4976-4982 (1992); 
Giovannangeli et al; Proc. NaU. Acad. Sci., 89: 8631-8635 (1992); Moser and Dervan! 
Science, 238: 645-650 (1987); McShan et al, J. Biol. Chem.. 267:5712-5721 (1992); 
Yoon et al, Proc. Natl. Acad. Sci., 89: 3840-3844 (1992); Blume et al. Nucleic Acids 
Research, 20: 1777-1784 (1992); Thuong and Helene, Angew. Chem. Int. Ed. Engl. 
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32: 666^690 (1993); Esciide et al, Proc. l>lati: Acad. Scf., 93:4365^369 (1996); and 
the like.' Conditions for annealing single-stranded or diipldc ^gs to their smgle- 
strarided or duplex complenients iarie well khowh, eJg. Ji et al; *Anali'Chem. 65:>1323- 
13^8 (1993)1 editor et aiiru;^ 
5 the advantage of hot i^'quirihg a'"strippihg" ieaction'with polymerase tb expose the ' ' 
tiag for annealing to its * ' ' • ' ^ ^ ' ■ * >^ n:iiVK.i ...r i < ^. 

^ Prtfefably/oligohikclebtide'tags of the invention employing triplex ^ ^ <^ - 
hybridization are double s^^hded DN A and 'the corresponding tag complementsare 
single sti^dedi ' More prefer^b ih place of cytosine in the 

10 tig 'com|ilements in order to brba^^ the range of pH stability of the triplex^fomied 
between a tajg arid its cbmpleniOTtV Preferred cond for forriung triplexes are ^ 
fiiliy disclosed in the above reference • > . 

concehtratecl salt'solutidn, e:g: 1 .0 M NaQ^l ;0 M potassium acetate; or the like, at r 
pH below 5.5 ( or 6.5 if 5-methylcyt6sine is employed). Hybridization temperature 

1 5 depends on theiength and Composition of the tag; however; for an 1 8-20rmer lag of ^ 
longer, hybridization at room temperature is adequate: Washes may-be conducted 
with less concentrated salt solmiohs - e.g. 10 ihM sodium acetate, 100 mM MgCJv,' pH 
5:8, at room temperature. Tags may be eluted from their tag complements by 
- incubation in a similar salt solution at pH 9.0. ' 

20 Minimally cross-hybridizing sets ofoligonucleotide tags that form triplexes 

may be generated by the computer program of Appendix Ic, or similiar programs. An 
exemplary set of double stranded 8-mer words are listed below in capital letters with 
the corresjx)nding complements in small letters. Each siich word differs from each of 
the other words in the set by three base pairs. ' * 

25 <^ ■ - . . . ^ . - 

Table V 

Exemplary Minimally Cross-Hvbridizinp 
Set of DoubleStranded 8-mer Tags 



C 1 

«/ 


-AAGGAGAG 


. 5' 


-AAAGGGGA 


5' 


-AGAGAAGA 


5' 


-AGGGGGGG 


3' 


-TTCCTCTC 


3' 


-TTTCCCCT 


3' 


-TCTCTTCT 


3' 


-TCCCCCCC 


3' 


-ttcctctc 


3' 


-tttcccct 


.3' 


-tctcttct 


. 3' 


-tccccccc 




-AAAAAAAA 


5' 


-AAGAGAGA 


- 5' 


-AGGAAAAG 


5' 


-GAAAGGAG 


3' 


1* I' I' 1^ 


3' 


-TTCTCTCT 


3' 


-TCCTTTTC 


3' 


-CTTTCCTC 


3' 


-ttt ctttt 


3' 


-ttctctct 


3' 


-tccttttc 


3' 


-ctttcctc 


i; ' 


-AAAAAGGG 


5' 


-AGAAGAGG 


5' 


-AGGAAGGA 


5' 


-GAAGAAGG 


3' 


-TTTTTCCC 


3' 


-TCTTCTCC 


3' 


-TCCTTCCT 


3' 


-CTTCTTCC 


3' 


-tttttCCC 


3' 


-tcttctcc 


3' 


-tccttcct 


3' 


-cttctrcc 


5' 


-AAAGGAAG 


5' 


-AGAAGGAA 


5' 


-AGGGGAAA 


5' 


-GAAGAGAA 


3' 


-TTTCCTTC 


3' 


-TCTTCCTT 


3' 


-TCCCCTTT 


3' 


-CTTCTCTT 


3' 


-rrtccttc 


3' 


-tcttcctt 


. 3' 


-tccccttt 


3' 


-Cttctctt 
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-V - . . Table yi 

Repertoire Size of Various Doubje Stranded Tags 
' That Form Triplexes wit h Their fag Complements ' 



1 



'i 1'. ! 



Oligonucleotid 
e 

Word" 
it.enflth 



Nucleotide'^. - ; ^v-:: *r,. l:'/^- 

Difference 

*' between^ • ' ^v: ---^M^irnal Sizfc- ^' 

. r-...., Oligonucleotides of MiniinaJly. .^^ v ;vSi?cof ... 

ofMimmaily Cross- * Rcpenoirc ' 

Hybridizing : - 



Crossr , I 



Hybridizing Set 



Set 



' iwilh Four 
Words 



Size of 
Repertoire with 



Five Words 



6 
8 
10 

15 

i 

20 
20 
20 



I I 



^ \ y 
3 

3 

5 

* ♦ 

8 

JO 



8 
8 
16 

92 
765 

92- 
22 



4096 



,1 1 



4096 
6.5 X. 10 
4096 



4 



3.2 X 10^ 
3.2; X I O'^ ' 
1.05 X 10^ 



Preferably, repenoires of double stranded oligonucleotide tags of the invention 
contain at least 10 members; more preferably, repertoires of such tags contain at least 
1 00 menibers. Preferably, words are between 4 and 8 nucleotides in length for 
combinatorially synthesized double stranded oligonucletide tags, and oligonucleotide 
tags are between i 2 and 60 base pairs in length. More preferably, such tags are 
between 1 8 and 40 base pairs in length. 

Solid Phase Supports 
Solid phase supports for use with the invention may have a wide variety of 
forms, including microparticles. beads, and membranes, slides, plates, micromachined 
chips, and the like. Likewise, solid phase supports of the invemion may comprise a 
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25 



:T''*^^..^^^^°^;?°'?'P°,?i^°"f.' incluciing glass, jjIm^^ • 

employed such that each has a unifpnn coating, or population, of complementarv 
5 sequences of thp same tag (and no other), or a single or a few supports are ernoloved 
with spatially diMj^te regidiB M^^^ uniform coating, or population, of 

comiiiementiMy>quencM todie.sanie tag (and no other). In the latter embodiment, 
the area of the f«!Si9ns may vary acx^^^ iiuajiy,']the 
regions range in aireafroin severai'mn\ e.g. 3-5, to several hundi^'um2" e g ' 100- 
• " , r^eferably, such regions are spau^lly discrete so that signals generated by 

events, e.g. fluorescent emissions, at adjacent regions can be resolved by the detection 
system being employed. In some applications, it may be desirable to have regions 
with unifomi coatings of more than one tag complement, e.g. for simultaneous 
sequence analysis, or for bringing separately tagged molecule's into close proximity. 
1 5 Tag complements may be used with the solid phase support that they are 

synthesized on, or they may be separately synthesized and attached to a solid phase 
support for use, e.g. as disclosed by Lund et al. Nucleic Acids Research, 1 6: 1 086 1 - 
10880 (1988); Albretseri et al. Anal. Biochem., 18,9: 40^50 (1990); W6lf et al. Nucleic 
- Acids Research: i5: 291 1-2926 (1987); or Ghosh 
5353-5372 (1987). Preferably, tag complements are sjoithesized on and used with the 
same solid phase support, which may comprise a variety oif forms and include a 
variety of linking moieties. Such supports may comprise microparticles or arrays, or 
matrices, of regions where uniform populations of tag complements are synthesized. 
A wide variety of microparticle supports may be used with the invention, including 
microparticles made of controlled pore glass (CPG), highly cross-linked polystyrene, 
acrylic copolymers, cellulose, nylon, dextran, latex, polyacrolein, and the like, 
disclosed in the following exemplary references: Meth. Enzymol., Section A, pages 
1 1-147, vol. 44 (Academic Press, New Yoric, 1976); U.S. patents 4,678,814; 
4,413,070; and 4,046;720; and Pon, Chapter 19. in Agrawal, editor. Methods in 
Molecular Biology, Vol. 20, (Humana Press, Totowa, NJ, 1993). Microparticle 
supports further include commercially available nucleoside-derivatized CPG and 
polystyrene beads (e.g. available from Applied Biosystems, Foster City. CA); 
derivatized magnetic beads; polystyrene grafted with polyethylene glycol (e.g.. 
TentaGelTM, Rapp Poiymere, Tubingen Germany); and the like. Selection of the 
support characteristics, such as material, porosity, size, shape, and the like, and the 
type of linking moiety employed depends on the conditions under which the tags are 
used. For example, in applications involving successive processing with enzymes, 
supports and linkers that minimize steric hindrance of the enzymes and that facilitate 
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access to substrat^ are prefeired. Other important factors to be considered in selecting 
the most appropriate microparticle support include size uniforaiity; efficiency as a 
synthesis support, degree to which surface area known, and optical properties, e.g. as 
explain more fiilly below, cl?M;smooA,beads provide.instrumentational advantages: .u i 
whOT hMcUing large numbers of h^ . . 

Exemplary linking n^oleties for at^ . , . 

microparticle surfeces, are disclosed jn,Pon et al, Biotechniques,. 6:768-775 (1988); S : . 
Webb, U.S. patent 4,659,774;Baraii>;,et.d,..totemational,p . i., ,; 

PCT/US91/p6103; Brown.et al, J. Chem. Soc, Commun.. 1989:;8.9J.-:893-Damhaet 
al, Nucleic Acids Research, J 8: 3813.382 1 (1990); Beattie et al. Clinical Gheratstry, 
39: 71 9-722 (1 993); Maskos and Southe^, Nucleic. Acids Research, .20: .1679-1,684 < 
(1992); and the, like. , - ,, j, - 

I ■ r 

As mentioned above, tag complementsj^^ also be synthesizedipn a single !' 
(ora few) solid phase supp9rt , to form an array of regions unifoimly coated withftag 
1 5 complements. That is, within each region in.such an arr?y,the same^tag complement 
is synthesized. Techniques for synthesizing such arrays are disclosed in McGall et al. 
International application PCT/US93/03767; Pease et al, Proc. Natl. Acad. Sci., 91 : 
5022-5026 ( 1 994); Southern and Maskos, International, application , ; 
.PCT/GB89/01 1 14; Maskos and Southern (cited above); Southern et al. Genomics, 13: 
20 , , 1 008- 1 0 1 7 ( 1 992); and Maskos and Southern, Nucleic Acids Research, 2 1 : 4663- 
4669(1993). 

Preferably, the invention is implemented with microparticles or beads 
uniformly coated with complements of the same tag sequence, Microparticle supports 
and methods of covalently or noricovalently linking oligonucleotides to their surfaces .,. 
25 are well kno\yn, as exemplified by the following references: Beaucage and Iyer (cited 
above); Gait, editor, Oligonucleotide^Synthesis: A Practical Approach (IRL Press, 
Oxford, 1984); and the.references cited above. Generally, the size and shape of a 
microparticle is not critical; however, microparticles in the size range of a few, e.g. 1 - 
2, to several hundred, e.g. 200-1000 jun diameter are preferable, as they facilitate the 

30 construction and manipulation oflarge repertoires ofoligonucleotide tags with 
minimal reagent and sample usage. 

In some preferred applications, commercially available controlled-pore glass 
(CPG) or polystyrene supports are employed as solid phase supports in the invention. 
Such supports come available with base-labile linkers and initial nucleosides attached, 
e g Applied Biosystems (Foster City, CA). Preferably, microparticles having pore 
size between 500 and 1000 angsu-oms are employed. 

In other preferred applications, non-porous microparticles are employed for 
their optical properties, which may be advantageously used when tracking large 
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numbers of microparticles oh planar supports, such as a microscope slide. ' * 
Panicularly preferred non-porous microparticles are the glycidal methacrylate (GMA) 
beads available from Bangs laboratories (Cime^^^^^ Siich microparticles are -^^^ 
useful in a variety of sizes and' derivatized with a variety of linkage'groiips'for 
synthesizing tags Pi^fei^bly" for riiassively parallel J . c' 

manipulations of tagged rriicropairticles, '5 ^m'diaitieter GMA beads are employed. 



a: 



tii,itrV,. I*,-' 



I * 



' Attkchihig Tags to Polvhucleotidef^ = i m - : r\ r> 



For Sorting onto Solid Phase Supports 



An important ibpecf of the invention is the sortirigk^ 
populations of polynucleotideis, e.g. from a cDN A library, to ihicrbparticles or to 
1 5 separate regions on a solid phase support such that each niicroparticlebr region has- 
substantially only one kind of polynucleotide attached. This objective is 
accomplished by insuring that substantially all different polynucleotides have 
different tags attached. This condition, in turn, is brought about by taking a sample of 
. the fiill ensemble of tag-polyniicleotide conjugates for analysis: (It is acceptable that 
20 identical polynucleotides have different tags, as it merely fesiilts in the same 

polynucleotide beiiig operated on or analyzed twice iri'two different location^'.) Such 
sampling can be carried out either overtly— for example, by taking a small volume 
from a larger mixture— after the: tags have been attached to the polynucleotides, it can 
be carried out inherently as a secondary effect of the techniques used to process the 
25 polynucleotides and tags, or sampling can be carried out both overtly and as an 
inherent part of processing steps. 

Preferably, in constructing a cDNA library where substantially all different 
cDNAs have different tags, a tag repertoire is employed whose complexity, or number 
of distinct tags, greatly exceeds the total number of mRNAs extracted from a cell or 
30 tissue sample. Preferably, the complexity of the tag repertoire is at least 1 0 times that 
of the polynucleotide population; and more preferably, the complexity of the tag 
repertoire is at least 100 times that of the polynucleotide population. Below, a 
protocol is disclosed for cDNA library construction using a primer mixture that 
contains a full repertoire of exemplary 9-word tags. Such a mixture of tag-containing 
35 primers has a complexity of 8^, or about 1 .34 x 10^. As indicated by Winslow et al. 
Nucleic Acids Research, 19: 3251-3253 (1991), mRNA for library construction can 
be extracted from as few as 10-100 mammalian cells. Since a single mammalian cell 
contains about 5x10^ copies of mRNA molecules of about 3.4 x 10"* different kinds. 
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by standard techniques one can isolate the mRNAirom about 100 cells, or : 
(theoretically) about 5 x^O', mRNA moleculej^, Comparing tf^ numter.to the , 
complexi^ of the primer mixture shows that without any additional steps, and even 
assuming that mRNAsi are,converted mtq cD^As with perfec^efTiqency (1% 
5 ffficienc:|r or IcM i&more.a^^ the cpNA libraiy cons^ructipnjprcitocol results in 
a population containing no more than,37% of the total numbq^pf different tags. .That 

of. *?.V»g repertoire, .TheprpbabUity, of, obtaining. a 
doiible imder these conditio^ is atout 5%, wWch is w^ the preferred range. .With 
10 «^RNA frpm 10 cells, the ftactipnpf the tag repertoire sampled is reduced;to only 
37%, even a^siminjg that;all processing steps take place at 100% efficiency. In 
,^^"^'^?*p!^f '^f processing steps for constructing cDNA .libranes are very 
* "™'^ of thur^jb^^ being that good library should contain about 10* cDNA clones 
from mRNA extracted from 10* mammalian cells. 
15 Use of larger amounts of mRNA in the above protocol, or for larger amounts 

of polynucleotides in general, where the number of such molecules exceeds the 
complexity of the tag repertoire, a tag-polynucleotide conjugate mixture potentially 
contains every possible pairing of tags and types of mRNA or polynucleotide. In such 
- cases, overt sampling may be implemented by removing a sample volume after a 
0 serial dilution of the starting mixture of tag-polynucleotide conjugates. The amount 
of dilution required depends on the amount of starting material' and the efficiencies of 
the processing steps, which are readily estimated. 

If mRNA were extracted from 10* cells (which would correspond to about 0.5 
Hg of poly(A): RNA), and if primers were present in about 10-100 fold concentration 
5 excess~as is called for in a typical protocol, e.g. Sambrook et al. Molecular Cloning, 
Second Edition, page 8.61 {10 ^L 1.8 kb mRNA at 1 mg/mL equals about 1.68 x 10'" 
; moles arid 10 #L18-mer primer at 1 mg/mL equals about 1 .68 x 1 0'' moles], then the 
total number of tag-polynucleotide conjugates in a cDNA library would simply be 
equal to or less than the starting number of mRNAs, or about 5 x 1 0" vectors 
0 containing tag-polynucleotide conjugates-again this assumes that each step in cDNA 
construction-first sttand synthesis, second su^d synthesis, ligation into a vector- 
occurs with perfect efficiency, which is a very conservative estimate. The actual 
number is significantly less. 

If a sample of n tag-polynucleotide conjugates are randomly drawn fix)m a 
5 reaction mixture~as could be effected by taking a sample volume, the probability of 

drawing conjugates having the same tag is described by the Poisson distribution 

-A. ' 
P(r)=e" (X)7r, where r is the number of conjugates having the same tag and X=np, 

where p is the probability of a given tag being selected. If n=l 0* and p=l/(l .34 x 
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10*), then ^.=.00746 and P(2)=2.76 x 10*^ Thus, a sample of one nMllipniniolecules 
gives rise to an expected number of doubles well within the preferred rangei^ Such a 
sample is -readily obtained as follows: Assume thatrthe,5.x 10 niRNASjarej.perfectly 
cpnwited 5 x lO" vectbre with tag-cDN A- conjugates as inserts arid that the 5 x 
5 ;10 vectors are in a reaction solution having a vbluine.of 100 \jX. Four l^rfpld serial 
dilutions may be carried out^fay transferring 10 ^1 from the original solution^into a; . 
vessel containing 90 ^l of an appropriate buffer, such as TE. This process may. be 
repeated'for three additional dilutions to obtain a 100 ^1 solution containing S x .10^ 
vector molecules per ^l. A 2 jil aliquot from this solution yields 10 vectors 

1 0 containing tag-cDNA conjugates as inserts. This sarhple is then amplified by straight 
forward transformation of a competent host cell followed by culturing. 

Of course, as mentioned above, no step in the above process proceeds with 
perfect efficiency. In particular, when vectors are employed to amplify a sample of 
tag-polynucleotide conjugates, the step of transforming a host is very inefficient. 

1 5 Usually, no more than 1 % of the vectors are taken up by the host and replicated. 
Thus, for such a method of amplification, even fewer dilutions would be required to 
obtain^sample of 10^ conjugates.' ' ^ • . ' : . . - , L 

A repertoire of oligonucleotide tags can be conjugated to a population of 
* polynucleotides in a number of ways, including direct enzymatic ligation, ^ 

20 arhplificatibn, e.g: via'PCR, lisinjg priniers contatining the tag sequences, and the like. 
The initial ligating step produces a very large population of tag-polynucleotide 
conjugates such that a single tag is generally attached to many different 
polynucleotides! However, as noted above, by taking a sufficiently small sarhple of 
the conjugates, the probability of obtaining "doubles," i.e. the same tag on two 

25 different polynucleotides, cari'be rriade negligible. Generally, the larger the sample 
the greater the probability of obtaining a double. Thus, a design trade-off exists 
between selecting a large sample of tag-polynucleotide conjugates-which, for 
example, ensures adequate coverage of a target polynucleotide in a shotgun 
sequencing operation or adequate representation of a rapidly changing mRNA pool, 

30 and selecting a small sample which ensures that a minimal number of doubles will be 
present. In most embodiments, the presence of doubles merely adds an additional 
source of noise or, in the case of sequencing, a minor complication in scanning and 
signal processing, as microparticles giving multiple fluorescent signals can simply be 
ignored. 

35 As used herein, the term "substantially all" in reference to attaching tags to 

molecules, especially polynucleotides, is meant to reflect the statistical nature of the 
sampling procedure employed to obtain a population of tag-molecule conjugates 
essentially free of doubles. The meaning of substantially all in terms of actual 
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percentages of tag-molecule conjugates depends on how the tags are being employed 
Preferably, for nucleic acid sequencing, substantially all means that at least eiehtv 
percent of the polynucleotides have unique tegs atteched. More preferably,,it means 
that at least ninety percent of the polynucleotides have unique tags atteched. Still 
5 more preferably, it rneans that at least ninety-five percent of the polynucleotides have 
unique tegs atteched. And, most preferably, it means that at least ninety-nine percent 
of the polynucleotides have unique tegs anached. 

Preferably, when the population of polynucleotides consists, of messeneer 
RNA (mRNA), oligonucleotides tags may be attached by reverse transcribine the 
1 0 mRNA with a set of primers preferably containing complements of tee sequences 
An exemplary set of such primers could have the following sequence (SEO ID NO* 

1): , . . , ;■ . . ■ " . 



15 



25 



30 



.5' -mRNA- (A)n -3' ' 

[ T j 1 9GG. [ W , W , W , CJ qAC CAGCTG ATC - 5 ' -bi ot i n 



t 

1 'It 



where "[W,W, W.CJg" represents the sequence of an oligonucleotide tag of nine 
- subimits of four nucleotides each and "[W.W,W,C]" represents the subunit sequences 
20 listed above, i.e. "W" represents T or A. The underlined sequences identify an 

optional restriction endonuclease site that can be used to release the polynucleotide 
from attechment to a solid phase support via the biotin, if one is employed. For the . 
above primer, the complement atteched to a microparticle could have the form: 



5'- [G, W,W,W] gTGG-linker-microparticle 

After reverse transcription, the mRNA is removed, e.g. by RNase H digestion, 
and the second strand of the cDNA is synthesized using, for example, a primer of the 
following form (SEQ ID NO: 2): 

5'-NRRGATCYNNN-3' 



where N is any one of A, T, G, or C; R is a purine-containing nucleotide, and Y is a 
pyrimidine-containing nucleotide. This particular primer creates a Bsi Yl restriction 
35 site in the resulting double stranded DNA which, together with the Sal I site, 

facilitetes cloning into a vector widi, for example. Bam HI and Xho I sites. After Est 
Yl and Sal 1 digestion, the exemplary conjugate would have the form: 
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5'-RCGACCA[C,W,W,W]9GG[Tli9- cDNA -NNNR 

GGT(G,W,W,W]9'Ct:[Arii9- rDNA -NNNYCTAG-5' 



/ 1 " 



The polynucleotide-tag conjugates may then be manipulated using standard molecular 
biology techniques. For example, the above conjugate— which is actually a mixture- 
may be inserted into conunercially available cloning vectors, e.g. Stratagene Cloning 
System (La JoUa, CA); transfected into a host, such as a conunercially available host 
bacteria; which is then cultured to increase the number of conjugates. The cloning 
vectors may then be isolated using standard techniques, e!g. Sambrook et al* 



1 0 Molecular Cloning, Second Edition (Cold Spring Harbor Laboratory, New York, 
1989). Alternatively, appropriate adaptors and piimers may be erhployed so that the 
conjugate population can be increiased By PCR: * 

. Preferably, when the ligase-based method of sequencing is employed, the Bsi 
Yl and Sal I digested fragments are cloned into a Bam HI-/Xho I-digested vector 
1 5 having the following single-copy restriction sites (SEQ ID NO: 3): ' ^ 



I • ; * :i 



20 



5 V-GA GGATG CCTTTAT GGATCC'AcTCGAG ATCCCAATCCA- 3 ' 
. . ; Fokl ' / *' BamHI ■ X'hol 

. I' - . . . , * * ' ' 



( t • 



This adds the Fok I site which will allow initiation of the sequencing process 
discussed more fiilly below. 

Tags can be conjugated to cDNAs of existing libraries by standard cloning 
methods. cDNAs are excised from their existing vector, isolated, and then ligated into 
25 a vector containing a repertoire of tags. Preferably, the tag-containing vector is 

linearized by cleaving vsdth two restriction enzymes so that the excised cDNAs can be 
. ligated in a predetermined orientation. The concentration of the linearized tag- 
containing vector is in substantial excess over that of the cDNA inserts so that 
ligation provides an inherent sampling of tags. 
30 A general method for exposing the single stranded tag after amplification 

involves digesting a target polynucleotide-containing conjugate with the 5*->3' 
exonuclease activity of T4 DNA polymerase, or a like enzyme. When used in the 
presence of a single deoxynucleoside triphosphate, such a polymerase will cleave 
nucleotides from 3' recessed ends present on the non-template strand of a double 
35 stranded fragment until a complement of the single deoxynucleoside triphosphate is 
reached on the template strand. When such a nucleotide is reached the 5'— ►3* 
digestion effectively ceases, as the polymerase's extension activity adds nucleotides at 
a higher rate than the excision activity removes nucleotides. Consequently, single 
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stranded tags constructed with three nucleotides are readily prepared for loading onto 
, solid phase supports., - 

The technique may also be used to preferentiaUy methy late, interior Fqk I sites 
P^!,a target polynucleot while leaving a single Fok Isite at the tenninus of the, 
poijmucleotide Fpk.I site is rendered single stranded 

3 u«n8^.P?:'y™^^^ with deoxycytidine.triphps stranded portion of 

the fragment is then methylated, after ^which the single stranded tenninus is filled in 
,^ PPlymerase in the presence of all four nucleoside triphosphates, thereby 
regenerating the Fok I site, Clearly, this procedure can be generalized to 
10 endonucleases other than Fok I. 

After the oligonucleotide tags are prepared for specific hybridization, e.g. by 
rendering them single stranded as described above, the polyimcleotides, are mixed 
with microparticles containing the complementaiy. sequences of the tags under 
conditions that favor the formation of perfectly matched duplexes between the taes 

15 and their compleinents. There is extensive guidance jn the literature for creating these 
conditions. Exemplary references providing such guidance include Wetmur, Critical 
Reviews in Biochemistry ai^d .Molecular Biology, 26: 227-25?,( 1 991 ); Sambrook et 
al. Molecular Cloiiing: A Lalwratoiy Manual,.2nd Edition (Cold Spring Harbor 
- Laboratory, New York, 1 989); and the like. Preferably, the hybridization conditions 

20 are sufficiently stringent so that only perfectly matched sequences form stable 

duplexes. Under such conditions the polynucleotides specifically hybridized through 
their tags may be ligated to the complementary sequences attached to the 
microparticles. Finally, the microparticles are washed to reniove polynucleotides with 
unligated and/or mismatched tags. . 

25 When CPG microparticles conventionally employed as s>Tithesis supports are 

used, the density of tag complements on the microparticle surface is typically greater 
than that necessary for some sequencing operations. That is, in sequencing 
approaches that require successive treatment of the attached polynucleotides with a 
variety of enzymes, densely spaced polynucleotides may tend to inhibit access of the 

30 relatively bulky enzymes to the polynucleotides. In such cases, the polynucleotides 
are preferably mixed with the microparticles so that tag complements are present in 
significant excess, e.g. from 10:1 to 100:1, or greater, over the polynucleotides. This 
ensures that the density of polynucleotides on the microparticle surface will not be so 
high as to inhibit enzyme access. Preferably, the average inter-polynucleotide spacing 

35 on the microparticle surface is on the order of 30-100 nm. Guidance in selecting 

ratios for standard CPG supports and Ballotini beads (a type of solid glass support) is 
found in Maskos and Southern, Nucleic Acids Research, 20: 1679-1684 (1992). 
Preferably, for sequencing applications, standard CPG beads of diameter in the range 
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25 



of 20-50 |im are loaded with about 10-^ polynucleotides, iarid GMA beads of diameter 
in the range of 5-10 \xm are'loaded witK a few teiis of th6usahdWpolynucleotides» 
e.g. 4 X 10^ to 6 X 10^. . >. r . . t 

In the preferred embodiment, tag complements are synthesized on 
5 microparticles combinatorially; thus, at the end of the synthesis, one obtains a 

complex mixture of microparticles from which a sample is taken for loading tagged 
polynucleotides. The size of the sample of niicroparticles will depend on several 
factors, including the size of the repertoire of tag complements, the nature of the 
apparatus^for used for observing loaded microparticles— e.g. its capacity, the tolerance 
1 0 for multiple copies of microparticles with the same tag complement (i.e. "bead 
doubles"), and the like: The following table provide guidance regarding 
microparticle sample size, microparticle diameter, and the approximate physical 
dimensions of a packed array of microparticles' of various diameters. 



Microparticle diameter 5 10 fun 20 ^un 40 inn ' 



i 



I • • 1 ■ 



' ' f 



Max. no. 

polyhucieotides loaded 
atlperlO^sq; ' ' ■ ! • 3x10^ J.26 xlO^; . ^5x10^ 

angstrom 

Approxl area of - • ' i r . , - 

monolayer of 10^ . , 

microparticles .45 x .45 cm 1 x I cm 2x2 cm 4x4 cm *' 



20 The probability that the sample of microparticles contains a given tag complement or 
is present in multiple copies is described by the Poisson distribution, as indicated in 
the following table. ' • 



Table Vll 
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Number of 



1 ^.-J-J 



Fraction of 
repertoire of tag 



micropanicles in . complements. : 
sample (as fraction present m 

ofrepertoire^izeX " ' '•••■^ > > ■ sample^ ' ' 
m l-e"*". 



Fraction of 
micropanicles in 

» > • ' V ! r: • . , >:-:rT i'Sampte, carrying 
, Fraction of . same tag 
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High Specificity Sorting and Panning 



The kinetics of sorting depends on the rate of hybridizat ion of oligonucleotide 
tags to their tag complements which, in tiirn, depends on the complexity of the tags in 
- the hybridization reaction, Thus^a trade off exists between sorting rate and tag . .. 
complexity,. such that an increase in sorting rate may be achieved; at the cost of , < 
reducing the complexity of the tags involved in the hybridi2ation reaction. As 
explained below, the effects pf this, trade off may be ameliorated by Tpanning." 

Specificity of the hybridizations may be increased by taking a sufficiently 
small sample so that both a high percentage of tags in the sample are unique and the 
nearest neighbors of substantially all the tags in a sample differ , by at least two words. 
This latter condition may be met by taking a sample that contains a number of tagr 
polynucleotide conjugates that is about 0.1 percent or less of the size of the repertoire 
being employed. For example, if tags are constructed with eight words selected from 
Table II, a repertoire of 8*, or about 1 .67 x 1 0^, tags and tag complements are 
produced. In a library of lag-cDNA conjugates as described above, a 0. 1 percent 
sample means that about 1 6,700 different tags are present. If this were loaded directi v 
onto a repertoire-equivalent of microparticles, or in this example a sample of 1 .67 x 
10*7 microparticles, then only a sparse subset of the sampled microparticles would be 
loaded. The density of loaded microparticles can be increase-for example, for more 
efficient sequencing«by undertaking a "panning" step in which the sampled tag- 
cDNA conjugates are used to separate loaded microparticles from unloaded 
microparticles. Thus, in the example above, even though a "0. 1 percent" sample 
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contains only 1 6,700 cDN As, ihe sampling arid panning steps may be repeatesd until 
as niahy loaded niicropaitid^^ * i ri/r : j*: 

A panning step may be implemented by providing a sample of tag-cDNA ' 
conjugates each of which contains a capture moiety at an end dppbsite, or distalno, 
5 the oligonucleotide tag. Preferably? the ca^^ typV wKicK 'can be ' 

released from tHe tag-cOl^A conjugiates/so tlia^ the lag-^DNA* conjugates can be 
sequenced with a single-Base sequencing ^ Siicii mbietibs'may 

biotin, digoxigenin, or likeligands, a triplex binding region, or t^^^ 
such a capture moiety corhpris^s a biotiri component. Biotih may be Attached to 

10 cDNA conjugates by a number of standard technique^! if apprbpriate adapters*' " 
containing PGR primer bihdinjg sites air^ attacHedlo tag-cDNA conjugates, bibiin may 
be attached by using a bidf mylated primer iii aii amplification after sampling. * ' 
Aitemiatively, if the tag-cDN A conjugates are inseit^'of cloning vectors, bibtin rnay be 
attached after excising the tkg-cDNA corijUgke by digestibh with an apprbpriate 

1 5 restriction enzyme followed by isoliatiori and filling in a prbtriidihg strand distal to the 
tags with a DNA polymerase in the presence of bibtinylated uridine triphosphate. 

After a tag-cDN A conjugate is captured,' it may be releatsed from the biotin 
moiet>' in a number of ways, siich as by a chemical linkage that 'is cleaved by " ' 
•reduction, e.g. Herman et al. Anal. Biochem., 156: 48-55 (1986), or that is cleaved 

20 photochemically, e:g. Olejhik et al, Niicleic Acids Research, 24: 361-366 (1996); or 
that is cleaved enzymatically by introducing a restrictibn site in the PCR primer. The 
latter embodiment can be exemplified by considering the library of tag-polynticleotide 
conjugates described above: ; • d , . i 

: ' » , . ' , - •* .... . - 

25 5'-RCGACCAiC,W,W,W]9GG[T]i9- cDNA -NNNR 

GGi'[G,W,'w,W]9CC[A]i9- rDNA -NNNYCTAG-5'' 

The following adapters may be ligated to the ends of these fragments to permit 
amplification by PCR: 

30 

5 ' - xxxxxxxxxxxxxxxxxxxx 

XXXXXXXXXXXXXXXXXXXXYGAT 

35 Right Adapter 



GATCZZACTAGTZZZZ2ZZZ2ZZZ-3 ' 
40 ZZTGATCA2ZZZZZZZZZZZ 
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, Left Adapter , : . 



ZZTGATCAZZ'ZZZZZ2Z2ZZ-'5 ' -biotin 



Left Primer 

wnere alt AOT" is a Spe I recognition site (which leaves a staggered cleavage 
ready for single base Mquehcing)! and the X's and'z nuclebtides selected so that 
the annealing and dissociation temperatures of the respective primers are 
approximately the same. After lipdonbf die adapters Md amplification hy PCR 
uising the biotinyiked primer,' 'the tag's of the conjugates are rendered single stranded 
by die exonuclease activity of 14 DNA polymerase Wd conjugate a^^^ with 
a sample of microparticles, e.g. a repertoire equivalent, with tag complements 
attached/ After annealing under stringent condhioiis (to minimize mis-^ttachmen^ 
tags), the conjugates are preferably ligated to their tag conipiemeiite 
micropiarticles are separated fi-om the unloaded micfop^cies by capture with 
avidinated majgnetic beads, or like capture technique. ' 

Returning to the example, this process r6suits in the acciimulation of about 
10,500 (=16,700 X .63) loaded microparticles with'diffcrent tags, which may be 
released from the magnetic beads by cleavage with Spe 1. By repeating this process 
40-50 times with new samples of microparticles and tag-cDNA conjugates^ 4-5 x 10^ 
cbNAs can be accumulated by pooling the released microparticles. ' The pooled ' 
microparticles may theii be simultaneously sequenced by a single-base sequencing 
technique: " ' ' ■ ■ . . 

Determining how many times to repeat thie sampling and panning steps--or 

riiore generally, determining how many cDNAs to analyze, depends on one's 
objective. If the objective is to monitor the changes in abundance of relatively 
common sequences, e.g. making up 5% or more of a population, then relatively small 
samples, i.e. a small fi-action of the total population size, may allow statistically 
significant estimates of relative abundances. On the other hand, if one seeks to 
monitor the abundances of rare sequences, e.g. making up 0.1 % or less of a 
population, then large samples are required. Generally, there is a direct relationship 
between sample size and the reliability of the estimates of relative abundances based 
on the sample. There is extensive guidance in the literature on determining 
appropriate sample sizes for making reliable statistical estimates, e.g. Roller et al. 
Nucleic Acids Research, 23:185-191 (1994); Good, Biometrika, 40: 16-264 (1953); 
Bunge et al, J. Am. Stat. Assoc., 88: 364-373 (1993); and the like. Preferably, for 
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monitoring changes in gene expression based on the analysis. of ^ series of cDNA 
libraries containing lO^-to 10* independent clones of 3.0-3.5 x 10^ different 
seqtiences,^ sample: of at least 10^ sequenceSjare accumulated for; wi^lysis of each^. 
library. More preferably, a sample of at least jip^, setju^ncesi are apcumulated for the ^ 
5 analysis of each library;: and? mpstpreferably, a;sample.of at least S x 10^ secjuences 
are accumulated for the analysis of each library. ; Alternj^tiyely, the number of 
sequences ^sampled is :preferably sufficient to estimate the. relative abundance of a 
sequence present atia frequency within the rang^,of;0.1% to.5,%.with a 95% , 
confidence limit ik> larger than 0:1;% of the popiulatipnsize. . 



rl^ :;ri:K::i: z\-i^r .vw Single Base E)N A Sequencing ^/v.-JH^i.iuv' ^ ^i: • 
; >' r The presem invention can :be employed with OTnv^^^^^ methods of DNA 
sequencings e g; as disclosed by Hultman et al, Nucleic Acids, Research, 1 7; 4937- 
4946 (1 989). However^ ifor parallel, or simultaneous, sequencing of multiple^ 

1 5 polynucleotides, a DNA sequencing; methodology is preferned that requires, neither; 
electrophoretic separation. of closely sized DNA fragments nor analysis of cleaved 
nucleotides by a separate analytical procedure, as. in peptide sequencing.^ Preferably, 
the methodology permits the stepwise 'identification of nucleotides, usually one.at a 
- lime, in a sequence ithroiigh successive cycles of treatment and detection. Such 

20 methodologies are referred to herein as "single base":sequencing methods.. Single, 
base approaches are disclosed in the following references: Cheeseman, U.S. patent 
5,302,509; Tsien et al. International application WO 91/06678; Rosenthal et al, ^ 
International application WO 93/21340; Canard etal. Gene, 148: 1-6 (1994); and 
Metzkeret al. Nucleic Acids Research, 22: 4259-4267 (1994). 

25 'A "single base" method of DNA sequencing which is suitable for use with the 

present invention and which requires no electrophoretic separation of DNA fragments 
is described in International application PCT/US95/03678. Briefly, the method 
comprises the following steps: (a) ligating a probe to an end of the polynucleotide 
having a protruding strand to form a ligated complex, the probe having a . 

30 complementary protruding strand to that of the polynucleotide and the probe having a 
nuclease recognition site; (b) removing imligated probe fi-om the .ligated complex; (c) 
identifying one or more nucleotides in the protruding strand of the polynucleotide by 
the identity of the ligated probe; (d) cleaving the ligated complex with a nuclease; and 
(e) repeating steps (a) through (d) imtil the nucleotide sequence of the polynucleotide. 

35 or a portion thereof, is determined. 

A single signal generating moiety, such as a single fluorescent dye, may be 
employed when sequencing several different target polynucleotides attached to 
different spatially addressable solid phase supports, such as fixed microparticles, in a 
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parallel sequencing operation. This may be accomplished by providing four sets of 
probes that are applied sequentially to ftcj. pliirality of target polynucleotides on the 
different microparticles. An exemplary set of such probes are shown below: 



•1 f/^, t 

Set I 



•-1 



Set 2 



Set 3 



Set 4 



ANNNN. . .NN 



dANNNN. .;.NN,/ , , . dAN^NR. ; wNN-^ , dANNNN. . . NN 

N. . .NNTT. . .T* d N.;,NNTT.;.t' ■ N. ..NNTT:..T 



N. . .NNTT. . .T 



dCNNNN. ..NN . . ^CNNNK.. .NN^, - dCNNNN . .^.NN,. . pdCNNNN : .,NN ' 

N. . .NNTT. . :t N. . .NNTT. . .T* N. . .NNTT.'. .T 



dGNNNN . . . NN 



N.". .NNTT. . .T 



N... .NNTT. . .T 



r.. , ^ ' f, 



■ . 1 ( 



dGNNNN NN 



, . GNNNN;';,NN.i 
N. . .NNTT.:.T^ 



— > • dGNNNN / It .;NN ';- '' 
N. . .NNTT. . .T* N...NNTT...T 



dTNNNN. . .NN 



N. . .NNTT: ;/.T ". f I- '"N. . ..NNTT. . .T 



dTNNNN . . .NN-'^' - 'TONNN. ' .NN' 



N. . .NNTT. . .t' 



'N. . .NNtT.\ .T* 



where each of the listed probes represents a mixture of 43=64 oligonucleotides such 
that the identity of the 3' terminal nucleotide of the top strand is fixed and the other 
positions in the.protruding strand are filled by every 3-mer permutation of nucleotides, 
or compjexity reducing analogs. The listed probes are also shown with a single 
stranded poly-.T, tail with a signal generating moiety attached to the tenninal thymidine, 
shown as ri*". The "d" on the unlabeled probes designates a ligation-blocking moiety 
or absense of 3'-hydroxyl, which prevents unlabeled probes from being ligated. 
Preferably, such 3'.tenninal nucleotides are dideoxynucleotides. In this embodiment, 
the probes of set lare first applied to the plurality of target polynucleotides and treated 
with a ligase so that target polynucleotides having a thymidine complementary to the 3' 
terminal adenosine of the labeled probes are ligated. The unlabeled probes are 
simultaneously applied to minimize inappropriate ligations. The locations of the target 
polynucleotides that form ligated complexes witii probes terminating in "A" are 
identified by the signal generated by the label carried on the probe. After washing and 
cleavage, the probes of set 2 are applied, hi this case, target polynucleotides forming 
ligated complexes with probes terminating in "C" are identified by location. Similarly, 
the probes of sets 3 and 4 are applied and locations of positive signals identified. This 
process of sequentially applying, the four sets of probes continues until the desired 
number of nucleotides are identified on the target polynucleotides. Clearly, one of 
ordinary skill could construct similar sets of probes that could have many variations, 
such as having protiiiding strands of different .lengths, different moieties to block 
ligation of unlabeled probes, different means for labeling probes, and the like. 
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A pparatus for Sequencing Populations of Polynucleotides 
An objective of the invention is to sort identical molecules, particularly 
polynucleotides, onto the sxirfaces of microparticles by the specific hybridization of 
tags and their complements. Once such sorting has taken place, the presence of the 
5 molecules or operations performed on them can be detected in a number of ways 
depending on the nature of the tagged molecule, whether microparticles are detected 
separately or in "batches,** whether repeated measurements are desired, and the like. 
Typically, the sorted molecules are exposed to ligands for binding, e.g. in drug 
development, or are subjected chemical or enzymatic processes, e.g. in polynucleotide 

1 0 sequencing. In both of these uses it is often desirable to simultaneously observe 

signals corresponding tO; such events or processes on large nimibers of microparticles. 
Microparticles carrying sorted molecules (referred to herein as "loaded" 
microparticles) lend themselves to such large scale parallel operations, e.g. as 
demonstrated by Lam et al (cited above). 

1 5 Preferably, whenever light*generating signals, e.g. chemiluminescent, 

fluorescent, or the like, are employed to detect events or processes, loaded 
microparticles are spread on a planar substrate, e.g. a glass slide, for examination with 
a scanning system, such as described in International patent applications 
* PCT/US9 1/092 17, PCT/NL90/00081, and PCT/US95/01886. The scanning system 

20 should be able to reproducibly scan the substrate and to define the positions of each 
microparticle in a predetermined region by way of a coordinate system. In 
polynucleotide sequencing applications, it is important that the positional 
identification of microparticles be repeatable in successive scan steps. 

Such scarming systems may be constructed from commercially available 

25 components, e.g. x-y translation table controlled by a digital computer used with a 
detection system comprising one or more photomultiplier tubes, or alternatively, a 
CCD array, and appropriate optics, e.g. for exciting, collecting, and sorting 
fluorescent signals. In some embodiments a confocal optical system may be 
desirable. An exemplary scanning system suitable for use in four-color sequencing is 

30 illustrated diagranunatically in Figure 5. Substrate 300, e.g. a microscope slide with 
fixed microparticles, is placed on x-y translation table 302, which is connected to and 
controlled by an appropriately programmed digital computer 304 which may be any of 
a variety of commercially available personal computers, e.g. 486-based machines or 
PowerPC model 7100 or 8100 available form Apple Computer (Cupertino, CA). 

35 Computer software for table translation and data collection functions can be provided 
by commercially available laboratory software, such as Lab Windows, available from 
National Instruments. 



35 



wo 97/13877 



rPCT/US96/16342 



Substrate 300 and table, 302 are pperationally associated with microscope 306 
having one or more objective lenses 308 which are capable of collecting and 
delivering light to microparticles fixed to substrate:300. Excitatioii.beam 3 1 0 from 
light source 312, which is preferably' a lasdr, is directed to beam splitt^^^^ a 

5 dichroic mirror, which re-directs the beam through microscope. 306 and objective lens 
308 which, in turn, focuses the beam onto substrate 300. Lens 308 collects 
fluorescence 316 emitted from the microparticles and directs it through beam splitter 
3 1 4 to signal distribution optics 31 8 which, in turn, directs fluorescence to one or 
more suitable opto-electronic deivices for conyertiiig some, fluorescence characteristic, 

0 e.g. intensity, lifetime, or the like, to an electrical signal. Signal distribution optics 
3 1 8 may comprise a variety of components standard in the art, such as bandpass 
filters, fiber optics, rotating mirrors, fixed position mirrors and lenses, diffraction 
gratings, and the like. As illustrated in Figure 2, signal distribution optics 3 1 8 directs 
fluorescence 316 to four separate photomultiplier tubes„33p, 332,.334, and 336, , 

5 Nvhose output is then directed to pre-amps and photon counters 350, 352,, 354, and 
3^6- The output of the photon counters is collected by computer 304, yyhere jt can be 
stored; analyzed, and viewed on video 360. Alternatively, signal distribution optics 
3 1 8 could be a diffraction grating which directs fluorescent signal 3 1 8 onto a CCD 
* array. 

The stability and reproducibility of the positional localization in scanning will 
detemiine, to a large extent, the resolution for separating closely spaced 
microparticles. Preferably, the scanning systems should be capable of resolving 
closely spaced microparticles, e.g. separated by a particle diameter or less. Thus for 
most applications, e.g. using CPG' microparticles, the scanning system should at least 
have the capability of resolving objects on the order of 10-100 ^m. Even higher 
resolution may be desirable in some embodiments, but vyith increase resolution, the 
time required to fully scan a substrate will increase; thus, in some embodiments a 
compromise may have to be made between speed and resolution. Increases in 
scanning time can be achieved by a system which only scans positions where 
microparticles are known to be located, e.g from an initial full scan. Preferably, 
microparticle size and scanning system resolution are selected to permit resolution of 
fluorescently labeled microparticles randomly disposed on a plane at a density 
between about ten thousand to one hundred thousand microparticles per.cm2. 

In sequencing applications, loaded microparticles can be fixed to the surface 
of a substrate in variety of ways. The fixation should be strong enough to allow the 
microparticles to undergo successive cycles of reagent exposure and washing without 
significant loss. When the substrate is glass, its surface may be derivatized with an 
alkylamino linker using commercially available reagents, e.g. Pierce Chemical, which 
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in tuni may cross-linkei^ to avidin, again using tbnventioni^I chemistries,^ 
ah ayidihated surface. Biotin moieties can 'te introduced to tKe loaded ihicrbparticles 
in a number of ways. For ex2firriple/a fraction; e:g. 10-15^1^ ' '^' ^ 

vectors used' to attach tags to pdlynuclebtides are engihedr^ t6' '(X)htain a unique 
5 restnction site (providing sticky endis on digestion) immediately adjacent to the 
polynucleotide insert aV an end of th[e polynucleotide opposite of the tag. The site is 
excised with the polynucleotide and tag for loading onto ihicrbparticles. "After 
1 oadihg, about 1 6- 1 5 percent the loaded polyriiicleotide^ will pos^eis the unique 
restriction site distzU from the micropai^^ surfice. After digestion wifti tiie ' '^' 
10 associated restrictibri endohuciease, an appfbpridte dbuble'sti^ded adaptor' ' 
containing a bibtin nioie^ is ligated to the sticky end. The resulting mictoparticles 
are then spread on the avidinated glass surface where they become fixed via the 

biotin-avidinlilikkges. ;^'"^^^^ '^■^«-.^^- • - ,.'^:-ir:.3r- 1 

Alternatively and preferably when sequencing by ligatidn is etnployed,1n the 
1 5 Initial ligation step a mixture of probes is applied to the loaded niicrbparticl6: a 

fraction of the probes contain a type'IIs restriction recogriitiori site, as required by the 
sequencing method, and k fraction of the probes have rib siich recbgiiition site, 'but 
inistead contain a biotin moiety at its non-ligating end. Preferably, the mixture 
* comprises about 10-15 percent of the biotinylated probe. ' . * ; . : 
20 In still another alternative, when DNA'-loaded micrbparticles are applied to a 

glass substrate, the DNA may nonspecifically adsorb to the glass surface upon several 
hours, e.g. 24 hours^ incubation to create a bond sufficiently strong tb permit repeated 
exposures to reagents and washes without significant loss of microparticles. ' 
Preferably, such a glass substrate is a flow cell, which niay comprise a channel etched 
25 in a glass slide. Preferably, such a channel is closed so that fluids may be pumped 
through it and has a depth sufficiently close to the diameter of the microparticles so 
that a monolayer of microparticles is trapped within a defined observation region. 

Identification of Novel Polynucleotides 
30 in cDNA Libraries 

Novel polynucleotides in a cDNA library can be identified by constructing a 
library of cDNA molecules attached to microparticles, as described above. A large 
fraction of the library, or even the entire library, can then be partially sequenced in 
parallel. After isolation of mRNA, and perhaps normalization of the population as 
35 taught by Soares et al, Proc. Natl. Acad. Sci., 91 : 9228-9232 (1994), or like 

references, the following primer may by hybridized to the polyA tails for first strand 
synthesis with a reverse transcriptase using conventional protocols (SEQ ID NO: 1 ): 
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5'-mRNA- [A]n -3' 

L i[.TJ 29>^lpl:imef • site) -GG [W; W/W, CJ'gACCAGCTGATC-S''^ 

where [ W,Wi, W,€]9 represent a tajg as described above, "AGCAGCfGATC" is kri ' 
5 optional sequence forming^ a restriction^ site in double straiitfied forniV aiid "primer site" 
is a sequence common to allimembers df the libfary thkt is latei- uied as a pnmer 
binding^site f6r;amplifying polynucleotides of interest by PCR. 

After reverse transcription ind second strand synthesis by con vention^ ^ 
techniques, ithe double stranded ifegments are inserted im 
0 described above arid' amplified: The amplified libraiyis then sampled arid the sample 
amplified. : The cloning vectors fi^oiri die amplified sample are isolated, and the tagged 
cDN A Augments excised andpurified; After rendering the teg single stranded with a 
polymerase as described above • the fi^gments are methylait^^ 
microparticles in accordance with the invention.' Preferably, as described above, the 
5 cloning vector is constructed 'so that the tagged cDNAs can be' excised with an 
endonuclease; such as Fok I, that will allow immediate sequencing by the preferred 
single base method after sorting and ligation to micropart 

Stepwise sequencing is then canied out simultaneously on the whole library, 
or one or more large fractions of the library, in accordance with the invention until a 
^sufficient number of nucleotides are identified on each cDN A for unique 
representation in the genome of the organism from which the library is derived. For 
example, if the library is derived fit)m mammalian mRNA then a randomly selected 
sequence 14-15 nucleotides long is expected to have unique representation among the 
2-3 thousand megabases of the typical mammalian genoriie. Of course identification 
of far fewer nucleotides would be sufficient for unique representation in a library 
derived from bacteria, or other lower organisms. Preferably, at least 20-30 ' 
nucleotides are identified to ensure unique representation and to permit construction 
of a suitable primer as described below. The tabulated sequences may then be 
compared to known sequences to identify unique cDN As. 

Unique cDNAs are then isolated by conventional techniques, e.g. constructing 
a probe fi-om the PCR amplicon produced with primers directed to the prime site and 
the portion of the cDNA whose sequence was determined. The probe may then be 
used to identify the cDNA in a library using a conventional screening protocol 

The above method for identifying new cDNAs may also be used to fingerprint 
mRNA populations, either in isolated measurements or in the context of a 
dynamically changing population. Partial sequence information is obtained 
simultaneously from a large sample, e.g. ten to a hundred thousand, or more, of 
cDNAs attached to separate microparticles as described in the above method. 
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^Example! 

An exemplary tag lil)ffflf^ is coiisBiic the chemically 

5 synthesized 9-\v6rff tags of nuclei 

3'-TGGC-[4(A,G,T)9]-CCCCp 

where "[^(A',G,T)9]" indicates a lagmixturew each tag consists of nine 4-mer 
words of A, G, and T; and "p" indicate a 5' phosphate. This mixture is ligated to the 
following right and left primer binding regions (SEQ ID NO: 4 and SEQ ID NO 5): 

-5'- AGTGGCTGGGCATCGGACCG 5'- GGGGCCCAGTCAGCGTCGAT 

TCACCGACCCGTAGCGp - - - u., ; GGGTCAGTCGCAGCTA 

LEFT RIGHT 

The right and left primer binding regions are ligated to the above tag mixture, after 
which the single stranded portion of the ligated structure is filled with DNA 
'polymerase then mixed with the right and left primers indicated below and amplified 
to give a tag library (SEQ ID NO: 6). 

Left PriiBttr 

25 

5 ' - AGTGGCTGGGCATCGGACCG 

5 ' - AGTGGCTGGGCATCGGACCG- [ ^ (A, G, T) 9) ^GGGGCCCAGTCAGCGTCGAT 
TCACCGACCCGTAGCCTGGC- (A, G, T) 9] -CCCCGGGTCAGTCGCAGCTA 

30 ~ 

CCCCGGGTCAGTCGCAGCTA- 5 ' 
Right Primer 

35 The underlined portion of the left primer binding region indicates a Rsr II recognition 
site. The left-most underlined region of the right primer binding region indicates 
recognition sites for Bsp 1201, Apa h and Eco O 1091, and a cleavage site for Hga I 
The right-most underlined region of the right primer binding region indicates the 
recognition site for Hga I. Optionally, the right or left primers may be synthesized 

40 with a biotin attached (using conventional reagents, e.g. available from Clontech 
Laboratories, Palo Alto, CA) to facilitate purification after amplification and/or 
cleavage. 
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primer binding site Ppu MI site 

(plasmid) -5 V -AAAAGGAGGAGGCCTTGATAGAGAGGACCT- 
, ... . ...-xTTTCCTCCTCCGGAA.CTATCTCTCCTGGAr 



primer binding site 

10 C -GTTTAAAC-GGAtCC-TCTTCCTCTTCCTCTTCC-3 * - (plasmid) 

-CAAATTTG-GCTAiSG-AGAAGGAGAAGGAGAAGG- 

: "t. t 

Bam HI site 

Pme 1 site 

15 

The plasmid is cleaved with Ppu MI and Pme I (to give a.Rsr Il-coinpatible eiid and a 
flush end so that the insert is oriented) and then methylated with DAM methylase. 
The tag-containing construct is cleaved with Rsr II and then ligated to the open 
plasmid. after which the conjugate is cleaved with Mbo I and Bam HI to permit 
20 ligation and closing of the plasmid. The plasmid is then amplified and isolated and 
- used in accordance with the invention. 

Example 3 

Changes in Gene Expression Profiles in Liver Tissue of Rats 

25 Exposed to Various Xenobiotic Agents 

In this experiment, to test the c^ability of the method of the invention to 
detect genes induced as a result of exposure to xenobiotic compounds, the gene 
expression profile of rat liver tissue is examined following administration of several 
compounds known to induce the expression of cytochrome P-450 isoenzymes. The 

30 results obtained fi-om the method of the invention are compared to results obtained 
from reverse transcriptase PCR measurements and inununochemical measurements of 
the cytochrome P-450 isoenzymes. Protocols and materials for the latter assays are 
described in Morris et al. Biochemical Pharmacology, 52: 781-792 (1996). 

Male Sprague-Dawley rats between the ages of 6 and 8 weeks and weighing 

35 200-300 g are used, and food and water are available to the animals ad lib. Test 
compounds are phenobarbital (PB), metyrapone (MET), dexamethasone (DEX), 
clofibrate (CLO), com oil (CO), and P*naphthofiavone (BNF), and are available from 
Sigma Chemical Co. (St. Louis, MO). Antibodies against specific P-450 enzymes are 
available from the following sources: rabbit anti-rat CYP3A1 from Human Biologies, 

40 Inc. (Phoenix, AZ); goat anti-rat CYP4A1 from Daiichi Pure Chemicals Co. (Tokyo, 
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Japan); monoclonal mouse anti-rat CYPl Al, monoclonal mouse anti-rat CYP2C1 1, 
goat anti-rat CYP2E1, and monoclonal mouse;anti-rat CYP2B1 from Oxford 
Biochemical Research^ tac^ XPxfon^ anti-rabbit IgG. 

rabbit anti-goat IgG and gpJit an^^mpiae JgQ) are.^avai Jackson 
5 ImmunoReseaB?h,Uboratpri?s,(West Gipye,'PA).\,., , .... . . 

Animals are administered either PB (100 mg/kg), BNF (100 mg/kg), MET 
(100 mg/kg), DEX (100 mg/kg), or CLO (250 mg/kg) for 4 consecutive days via 
intraperitoneal injection following a dosing regimen similar to that described by 
Wang et al. Arch. Biochem. Biophys. 290: 355-361 (1991). Animals treated with 
1 0 H2O and CO are used asxoritrpls. .Two hoiiirs following the last injection. (day 4), 
animals are killed, and the livers%e;rem6vied. Livers, are imrriediaiely frozen ancj 
stored at -70OC. 

Total RN A is prepared from frozen liver tissue using a modification of the 
method described by Xie et al, BiotecHniques, 1 1 : 326-327 (1991 ). Approximately 
1 5 1 00-200 mg of liver tissue is homogenized in the RN A extraction buffer described by 
Xie et al to isolate total kN A. The resulting RNA is reconstituted in 
diethylpyrocarbonate-treated water, quantified spectrophotometrically at 260 nm, and 
adjusted to a concentiation of 100 fig/ml. Total RNA is stored in 
- diethylpyrocarbonate-treated water for.up to 1 year at -70OC without any apparent 
20 degradation! Rt-PCR and sequencing are performed on samples from these 
preparations. 

For sequencing, samples of RNA corresponding to about 0.5 ^g of poly(A)+ 
RNA are used to construct libraries of tag-cDNA conjugates following the protocol 
described in the section entitied "Attaching Tags to Polynucleotides for Sorting onto 
Solid Phase Supports," with the following exception: tiie tag repertoire is constructed 
from six 4-nucleotide words from Table II Thus, the complexity of the repertoire is 
86 or about 2.6 x 1 05. For each tag-cDNA conjugate library constructed, ten samples 
of about ten thousand clones are taken for amplification and sorting. Each of the 
amplified samples is separately applied to a fixed monolayer of about 10^ 10 pm 
diameter GMA beads containing tag complements. That is, the "sample" of tag 
complements in the GMA bead population on each monolayer is about four fold the 
total size of the repertoire, thus ensuring tiiere is a high probability tiiat each of Uie 
sampled tag-cDNA conjugates will find its tag complement on the monolayer. After 
the oligonucleotide tags of the amplified samples are rendered single stranded as 
described above, tiie tag-cDNA conjugates of tiie samples are separately applied to the 
monolayers under conditions that permit specific hybridization only between 
oligonucleotide tags and tag complements forming perfectiy matched duplexes. 
Concentrations of the amplified samples and hybridization times are selected to 
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permit the loading 6f about 5 x l(H-to 2 x 10^ tag-cDNA conjugates on each bead 
where perfect matches occur: After ligation, 9-12 nucleotide portions of the attached 
cDNAs are determined in parallel by thb single base sequencing technique described 
by Brenner in International patent application PCT/US95/03678. Frequency 
5 distributions for the gene expression^profiles are assembled from the sequence 
information obtained from each of the ten samples. 

RT-PCRs of selected mRNAs corresponding to cytochrome P-450 genes and 
the constitutively expressed cyclophilin gene are carried out as described in Morris et 
al (cited above). Briefly, a 20 |liL reaction mixture is prepared containing Ix reverse 
0 transcriptase buffer (Gibco BRL), 10 nM dithiothreitol, 0.5 nM dNTPs, 2.5 oligo 
d(T)|5 primer, 40 units RNasin (Promega, Madison,: WI), 200 units RNase H-reverse 
transcriptase (Gibco BRL), and 400 ng of total RNA (in diethylpyrocarbonate-treated 
water). The reaction is incubated for 1 hour at 37^0 followed by inaclivation of the 
enzyme at 95^C for 5 miri. The resulting cDNA is stored at -20^C until used. For 
5 PCR amplification of cDNA, a 10 pL reaction mixture is prepared containing lOx 
polymerase reaction buffer. 2 mM MgCl2, 1 unit Taq DNA polymerase (Perkin- 
Elmer, Norwalk, CT), 20 ng cDNA, and 200 nM concentration of the 5* and 3' 
specific PCR primers of the sequences described in Morris et al (cited above). PCRs 
-are carried out in a Perkin-Elmer 9600 thermal cycler for 23 cycles using melting, 
0 annealing, and extension conditions of 940C for 30 sec, 56^C for 1 min., and 72^C 
for 1 min., respectively. Amplified cDNA products are separated by PAGE using 5% 
native gels. Bands are detected by staining with ethidium bromide. 

Western blots of Ae liver proteins are carried out using standard protocols 
after separation by SDS-PAGE. Briefly, proteins are separated on 1 0% SDS-PAGE 
5 gels under reducing conditions and immunoblotted for detection of P-450 isoenzymes 
using a modification of the methods described in Harris et al, Proc. Natl. Acad. Sci., 
88: 1407-1410 (1991). Protein are loaded at 50 ^g/lane and resolved under constant 
current (250 V) for approximately 4 hours at 2^C\ Proteins are transferred to 
nitrocellulose membranes (Bio-Rad, Hercules, CA) in 15 mM Tris buffer containing 
120 mM glycine and 20% (v/v) methanol. The nitrocellulose membranes are blocked 
with 2.5% BSA and inmiunoblotted for P-450 isoenzymes using primary monoclonal 
and polyclonal antibodies and secondary alkaline phosphatase conjugated anti-IgG. 
Immunoblots are developed with the Bio-Rad alkaline phosphatase substrate kit. 

The three types of measurements of P-450 isoenzyme induction showed 
substantial agreement. 
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APPENDIX la 

^ ^ , Exemplar y computer program for generating 
minimally cross hybridizing sets 
(single stranded tag/single stranded tag complement) 

Program minxh 

c - ■ ' ' ' ■' " ' • 

c 

integer*2 subl (6) , mset 1 ( 1000, 6) , mset2 ( 1000, 6) 

dimension nbase(6) 

c ... ... . . . , , . . 

c ■ * . . , ^ . , • 

write (*,*) 'ENTER SUBUNIT LENGTH' 
read (*, 100) nsub 
100 format (il) 

open {1, f ile= ' sub4 . dat * , f onn= ' formatted • , status^ ' new ' ) 

c 
c 

nset=0 

do 7000 ml=l, 3 
do 7000 m2=l,3 
do 7000 raS^l, 3 
do 7000 m4=l, 3 
subl (l)=ml 
subl (2)=m2 
subl {3)=m3 
subl (4)=m4 

c 
c 

ndiff=3 

c 
c 

c Generate set of subunits differing from 

c subl by at least ndiff nucleotides, 

c Save in msetl . 

c 

c 

do 900 j=l,nsub 
900 msetia, j)=subl(j) 

c 
c 

do 1000 kl«l, 3 
do 1000 k2=l,3 
do 1000 k3=l,3 
do 1000 k4=l,3 

c 
c 

nbase(l)=kl 
nbase(2)=k2 
nbase(3)=k3 
nbase (4 ) =k4 
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1200 

C 

c 



1 

3 



.1 If 



do 1200 j = l,nsiib' 
c '.c r.irf ,<suhl (;j};i}eqi 1 i.and..-. qbas^^( j ),r^nje 1 .or. 

subl ( j ) . eq. 2 , . aind. nbase ( j ) .rieT2 .or . 
su!blVjV>e'q:K3 ? ;WndV^TibiaseH[^ ) .ne. 3) then 
r^n^n+;i: ; - ^ r:- j.-i, ^ ^-r>.; ijC-i^r 
endif 
continue 



if (n,.ge.. ndif.f-) then 



c 
c 
c 
c 
c 
c 
c 



1100 

c 
c 

1000 

c 

c 



1325 



If number of mismatches 
is greater than or equal 
to ndiff then record 
subunit in matrix mset 



do 1100 i=l,nsub 

mset 1 ( j j , i ) =nbase ( i ) 
endif 



continue 



do 1325 j2=l,nsub ' 
mset2 (1, j2) =msetl (1, j2) 
mset 2 (2, j2 ) =msetl (2, j2) 



c 
c 
c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

c 
c 

1700 



npass=0 



continue 

)ck=npass+2 

npass=npass+l 



Compare subunit 2 from 
msetl with each successive 
subunit in msetl, i.e. 3, 
4,5, ... etc. Save those 
with mismatches .ge. ndiff 
in matrix mset2 starting at 
position 2. 

Next transfer contents 
of mset2 into msetl and 
start 

comparisons again this time 
starting with subunit 3. 
Continue until all subunits 
undergo the comparisons. 
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2 
2 

1600 

1625 

1500 

c 
c 
c 
c 

V-. 

c 
c 

2000 



c 
c 

7009 

700B 
7010 

120 
7000 

c 
c 

c 
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do 1500 in==npass+2, j j - . - 

n=0 ... 

do 1600L,j«l,nsub 

if (msetl (hpass+1, j ) .eq. 1 .and.msetl {m, j ) . ne, 1 . or . 
misetl {npa^s+1, j ) . eq. 2 . and.msetl (m, j ) . ne. 2 . or . 
msetl (npass+1, j) . eq. 3. and.msetl {m, j ) .ne.3) then 
n=n+ l ' . 

endif 

continue 

i f { n . ,ge .^ndi f f ) then 

: r.r.kk=kk+lK ■ 

>^:do 162& i=l,nsub 

' mset2fkk,i)=msetl{m,i) 
endif • 

cont'inue „ 

■ ' i,' '■ ' ' '. . • • ■ . 

i'r.-.^. ; ; ' kk is the number of subunits 

V. ^ ■ , stored in mset2 

4 ♦ * - . I 

Transfer contents of mset2 
y into msetl for next pass. 

do .2000 k=l, kk 

do ,2000 m=l,nsub . 

mset 1 (k,-m)=mset2 (k,mr 
if (kk. it . j j ) then • . • 

goto 1700 
endif ' 



nset=nset+l 
write (1, 7009) 

format { / ) 
do 7008 k=l, kk 

write (1, 7010) (mset 1 ( k, m) , m-1 , nsub) 
format.(4il) 
write(*,*) 

write{*, 120), kk,nset 

format ( Ix, 'Subunits in set= * , i5, 2x, ' Set No=M5) 
continue 
close (1) 



end 



46 



wo 97/13877 



PCT/US96/16342 



APPENDIX lb 

Exemplary computer program for generating 
minimally cross hybridizing sets 
(single stranded tag/single stranded tag complement) 



Program tagN 

c ' •^ •* : ' ' ■: 

c 

c Program. tagN generates minimally cross-hybridizing 

c* V sets of "suburii'ts given i) N — subunit length, and ii) 

c .an initial subunit sequence. ' tagN - assumes that only 

c 3 of the four natural nucleotides are used in the tags 

character*! subl (20) 

integer*2 mset (10000, 20) , nbase(20) 

c 

write (^;*) 'ENTER SUBUNIT LENGTH * 
;read'(*/lO0)nsub 
100 . format (i2) -'-^ ' 
c ...... 

c 

^: write:{*, *) 'ENTER SUBUNIT SEQUENCE* 
.; .1 read( * , 1 1 0 ) *( subl ( k)- k= 1 , nsub ) 
110 * . forrnat;(20al) 

c 
c 



c 
c 



ndiff=10 



c Let a=l c=2 g=3 & t=4 



do' 800 kk=l, nsub ' 
if (subl (kkj .eq. 'a' ) then 
mset(l,kk)=l . 

•endif " ; 

if!(subl (kk) .eq. •c.':) then 
mset(l>kk)=2 
endif 

if (subl (kk) .eq, 'g' ) then 
mset(l,kk)=3 
endif 

if (subl (kk) .eq. ' t • } then 
mset (1, kk) =4 
endif 

800 continue 

c 
c 

c Generate set of subunits differing from 

c subl by at least ndiff nucleotides. 



do 1000 ki=l,3 
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do 1000 k2=l,3 
do 1000 k3=l,3 

do. 1000 k4=l,3 
^ do 1000 k5=l,3 

do 1000 k6=l,3 
, . do 1000 k7=l,3 

do 1000 k8=l,3 
do 1000 k9=l,3 
do 1000 klO=l,3 

do 1000 kll-1,3 
• i ' ^. ' do 1000.:kl2'=l, 3 

do 1000 kl3=l,3 
do 1000 kl4=l,3 
• . .... ^ ^do' 1000^ klS^l, 3 
■ ' ; do iOOO *kl6=l, 3 

; . . ; , ': ,.^:^-/.dc, ipOO . kl7=i, 3 
~. * .*■■- . do' 1000 kl8=l,3 

" ' do 1000 kl9«l, 3 

do 1000 k20=l,3 

c 



nbase { 1 ) 


=kl 


nbase (2) 


= k2 


nbase (3) 


=k3 


nbase ( 4 ) 


=k4 


nbase ( 5 ) 


= k5 


nbase (6) 


=k6 


nbase (7 ) 


-k7 


nbase (8) 


=k8 


nbase (9) 


-k9 


nbase { 10 


), = klO 


nbase (11 


)«kll 


nbase { 12 


)=kl2 


nbase {13 


)=kl3 


nbase {14 


)=kl4 


nbase (15 


)=kl5 


nbase (16 


)=kl6 


nbase ( 17 


)=kl7 


nbase (18 


)=^kl8 


nbase (19 


)=kl9 


nbase (20 


)=k20 



do 1250 nn=l, j j 
n=0 

do 1200 j=l,nsub 

if (mset (nn, j) .eq.l .and. nbase { j ). ne . 1 .or. 
mset (nn, j ) .eq.2 .and. nbase ( j ) . ne . 2 .or. 
mset (nno ) .eq.3 .and. nbase ( j ) . ne. 3 .or. 
mset (nn, j ) .eq.4 .and. nbase ( j ) . ne . 4 ) then 
n=n+l 
endif 
continue 



if (n.lt.ndiff ) then 

goto 1000 

endif 
continue 



write ( *, 130) (nbase (i) , i=l, nsub) , j j 
do 1100 i=l,nsub 



1 
2 
3 



1200 

c 

c 



1250 
c 
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. . ' mset ( jj, i)=nbase (i) 
1100 continue 

w 

r . ! ' ' . ■ • / . ' - ■ f • ■ 

c* ■ - - ; ■ - ; . . ~- 

1000 continue 

c 

write(*, *) 
130 format (lOx, 20 (Ix, il),5x,i5) 

write ( * / * ) 

write (*, 120) j j 
120 ' format*(lx, 'Number of worcis=*,i5) 
c 

■end • . : ^ , . ^ 
c -^*********;***.************* 
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APPENDIX Ic 

Exemplary computer proeram for generating ,,. . .^of.;,? 

* 1 

minimall y cross hyhndizinp sets 
,;,.<dpuble,s^Mde4:tag/single stranded tag complement) 



Program 3tagN 



"1 '>fw*::;;^T 



Ml 



■ f Til * 



c 
c 
c 
c 
c 



c 
c 



100 

c 

c 



110 
c 

* 

c 
c 



800 



c '1-. •) 



' Program! '3tagN'*generate^ iriinimally ^cros^Kybrib^ ' ' ' " 
. sets p^of^^ duplex subunits .giyen..i) |.N-rsubunit .^^^^^ 
and ii)' an'ihitial homopurine sequence. 



character*l::subl'(-20)'' • • 
integer*2 msec (10000, 20) , nbase(20) 



i \ 



write (*/*)' 'ENTER SUBUNIT LENGTH' 
,,read(*, 100)nsubo 
format (i2) , 



wfite(*,*) -ENTER 'SUBUNIT SEQUENCE "a & g only' 

read(*, 110)..(subl.(k).,.k'=L,nsub) r ' 

format (20al ) 



* I . ^ ! " * i 



ndif f=10 



.Let a^^l and g-2 



" do 80O 'kk=l,nsub • 
if (subl (kk) .eq. 'a' ) then 
" m'set (1, kk)=l ' • 
endif . 

if (subl (kk) .eq. 'g' )' then 
mset (1, kk)«=2 
endif 

continue 



jj = l 

do 1000 kl=l, 3 
do 1000 k2=l,3 
do 1000 k3=l, 3 
. do 1000 k4=l,3 
do 1000 k5=l, 3 
do 1000 k6=l,3 
do 1000 k7=l, 3 
do 1000 k8=l,3 
do 1000 k9=l,3 
do 1000 kl0=l,3 

do 1000 kll=l,3 
do 1000 kl2=l,3 
do 1000 kl3=l, 3 
do 1000 kl4'=l,3 
do 1000 kl5=l,3 
do 1000 kl6=l,3 
do 1000 kl';=l,3 
do 1000 kl8=l, 3 
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do 1000 kl9=l,3 
do 1000 k20=l,3 



I 'J 



J- * V. 



} • 'J 



C 

c 



1 

2 
3 



'1200 
c 



1250 
c 



1100 
c 

1000 
c 

130 



120 

c 

c 



nbase. 
hbase 
nbase 
nbase 
nbase 
nbase 
' hbaise 
nbase 
nbase 
nbase 
nbase 
•hbaise 
nbase 
nbase 
nbase 
nbase 
' nbase 
nbase 
' ribase 
, nbas^, 



2-)=k2 

3) =k3 

4) =k4 

5) =k5 

6) =k6 

8) =k8 

9) =k9 

10) =klO 
ID^kll ^ 
l'2^=kl2^^ 

13) =kl3 

14) «kl4 

15) =kl5 

16) =kl6 

17) «kl7r 

18) =kl8 

l9")=k'r9'^^ ' 
2Q)rrk20.,. 



t ; t 



do •1250 nn^lv jj, 



I 1 



n=o ■ ■ • ... 
do 1200 j=l,nsub 

if (mset (nn, j ) .eq. 1 .and. nbase (j ). ne . 1 .or. 
, ' . .. mset (nn, jl v,eq,.2. .and. .nbase{3l .or. . 

mset (nn, j ) . eq . 3 .and. nbase ( j ) . ne . 3 .or. 

mset {nn, j 1 . eq . 4 .and. nbase ( j ) . ne . 4 ) then 

n-n+1 

endif 

... continue.,. ... ... . / . ^ 



if (n , It . ndif f ) then 

goto 1000 

endif 
continue 



J3 = 



write (' , 130) { nbase (i ) , i=l, nsub) , j j 
do 1100 i=l,nsub 

mset ( j j , i ) =nbase (i ) 
continue 

continue 
write ( * , * ) 

format (lOx, 20 ( Ix, il) ,'5x, i5) ■ 
wr.ite ( * , * ) 
write(*,120) jj 

f ormat ( Ix, • Number of words=',i5) 



end 
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SEQUENCE LISTING 

(1) GENERAL;:INFORMATION;-. I ^ ■. j- . i 1^ t^^t' ; rU u-/- 

(i) APPLICANT: David W. Martin, Jr. " 

(ii) TITLE-.OF-:INVENTION: : Measurement of 'Gene 'Express idtf-dro'f lies in 
Toxicity Determination . 

(iii) NUMBER Or'SEQUENCEsV "9 , -.'v.-^ . F'i V. I ; .r ; '.|.; 

{iyr^-CORRESPONDENCE .ADDRESS: v ... 

(A) ADDRESSEE: Stephen c/'Macevicz/'Lynx* Wrapeutics/ Inc 
-(B) STREET: 3.832 'Bay , Center , Place 

(C) CITY: Hayward" 

(D^-5TATE,r <Caiiiornia . ' vt^u'- rr ' -.. m -0 ■ j M\A 

(E) COUNTRY: USA ' ' 

-r (-F.)' -ZIP:-. i94545 ' :.>^'.-'vr- -fcvs i^'-f , ."^^; 

(V) COMPUTER READABLE FORM: . 

* (A)' MEDIUM':TYPE':'.3;5 'inch diskette'''^' ^ ' ■ ■•' »^ 
. (B) COMPUTER: IBM compatible 
'(C) OPERATING SYSTEM: Wi'ndoWs 3 . l"' ' " ' "* 

(D) SOFTWARE: Microsoft Word 5.1 



(vi) CURRENT APPLICATION DATA: 

(A) APPLICATION NUMBER:- - , ■ ' .i 

(B) FILING DATE: 

• (C) , CLASSIFICATION: ^' 

(vii) PRIOR APPLICATION DATA: 

(A) APPLICATION NUMBER: PCt/US96/09513 

(B) FILING' DATE: 06-JUN-96 ' " 



(vii) PRIOR APPLICATION DATA: 

(A) APPLICATION NUMBER: ' PCT/US95/1279i 
. (B) FILING DATE: 12-OCT-95 

(viii) ATTORNEY /AGENT INFORMATION: 

(A) NAME: Stephen C. Macevicz 

(B) REGISTRATION -NUMBER: 30, 285' 

(C) REFERENCE/DOCKET NUMBER: 813wo 

••■*(■' ■,••' 

(ix) TELECOMMUNICATION INFORMATION: 

(A) TELEPHONE: (510) 670-9365 
' (B) TELEFAX: (510) 670-9302 



(2) INFORMATION FOR SEQ ID NO: 1: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 11 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 
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1* 



-{x'i)- 'SEQUENCE DESCRIPTION: SEQ ID NO: 1: 



> CTAGTCGAGC-A. ..,: , . , 



^ (2) INFORMATION FOR SEQ ID NO: 2: 

>M i;y- SEQUENCE . CHARACTERISTICS:' >f;.^••i)• . . " o-. ' i.,.v .-ri - 
. (A) LENGTH: 11 nucleotides 
' ' ' (B) TyPEi*"-nuclei'c'-'acid ' /\ ^ '^u y • J'., ■'>■''.. 
(C; STRANDEDNESS: , ^single 
' ' (D) 'TOPOLOGY': ''linear^ i-:. - ':■.»..■•■} 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 2: 



• 1- • 'j-i^i 



. NRRGATCYNR.N . , . . 1 



1 



[2] INFORMATION* FOR SEQ 'lb NO: 3: 



(i) SEQUENCE .CHARACTERISTICS,: , 

(A) LENGTH: 38 nucleotides* 
- ' (BJ ^TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 



(xi) .SEQUENCE DESCRIPTION: SEQ ID NO: 3: 

' ■ ' .* ■■ . . ■ • , * , 

GAGGATGCCT TTATGGATCC ACTCGAGATC CCAATCCA 38 



{2r INFORM^i.TION FOR SEQ ID NO: 4: 

. (i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 20 nucleotides 

(B) -TYPE:: nucleic acid . 

(C) STRANDEDNESS: double 
'(D) TOPOLOGY; linear 



(xi;- SEQUENCE DESCRIPTION: SEQ ID NO: 4: 



AGTGGCTGGG CATCGGACCG 2C 



• (2) INFORMATION FOR SEQ ID NO: 5: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH; 20 nucleotides 

(B) TYPE: nucleic acid 
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(C) STRANDEDNESS: double 

(D.) . TOPOLOGY.:, -linear \ v. -- ^. 

Ui)'' SEQUENCE- DESCRI PTION:' SEQ ID* NO: 5': ' 



GGGGCCCAGT CAGCGTCGAT 

(2.) I^NFORMATION.rFOR/SEQ. ID N0: -6:.'i ^ IV^. - . 

(i) SEQUENCE CHARACTERISTICS: 

(AJ LENGTH: 20 nucleotides 

(B) TYPE:, nucleic acid 

(O- STRANDEDNESS : single" ^ " ^ ' 

(D) .TOPOLOGY; .linear . . 

(xi) 'SEQUENCE DESCRI PTIOnV'SEQ ID NO: 6: " 



ATCGACGCTG ' ACTGGGCCCC ' 



(2) INFORMATION FOR SEQ ID, NO.:- 7: 

(i)* SEQUENCE CHARACTERISTICS: ' ' " 
' . - , (A);.LENGm:::,62 nucleotides 
(B) TYPE: nucleic acid 
(G) ,STRANDEDNESS: double ' - 

(D) TOPOLOGY: linear 

(xi) SEQUENCE DESCRI PTION : -^SEQ ID NO: 7:' 

AAAAGGAGGA GGCCTTGATA GAGAGGACCT' GTTTAAACGG ^ATCCTCTTCC 
TCTTCCTCTT CC ' ' ^^'^ — . • * 
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I claim: 



1 . A method of determining the toxicity of a compound, the method comprising 
the steps of: ... * '^^'^^^ 

administering the compound to'a test prgam ... 

extracting a population of mRNA molecules from each of one or more tissues 



of the test organism; 

forming a separate population of cDNA mblecules from ea^^ population of 
mRNA molecules from the one or more tissues such that each cDN A molecule of a 
1 0 separate population has an oligonucleotide tag attached,'the oligonucleotide tags 
being selected from the same ininimally cross-hybridizing set; - - . 

separately sampling each population of cDNA niolecules such that; 
substantially all different cDNA molecules; within a separate po|)iJ^^^ have different 
oligonucleotide tags attached; . : ~ ; ^ su- ^r: ' ; ■ 

1 5 sorting the cDNA molecules of each sepiarate population by specifically 

hybridizing the oligonucleotide tags with their respective complements, the respective 
complements being attached as uniform populations of substantially identical 
complements in spatially discrete regions on one or more solid phase supports; 

determining the nucleotide sequence of a portion of each.of the sorted cDNA 
20 molecules of each separate population to forni a frequency distribution of expressed 
genes for each of the one or more tissues; and 

correlating the frequency distribution of expressed genes in each of the one or 
more tissues with the toxicity of the compound. - 

25 2. The niethod of claim 1 wherein said oligonucleotide lag and said complement 
of said oligonucleotide tag are single stranded. 

3. The method of claim 2 wherein said oligonucleotide tag consists of a plurality 
of subunits, each subunit consisting of an oligonucleotide of 3 to 9 nucleotides in 

30 length and each subunit being selected from the same minimally cross-hybridizing set. 

4. The method of clairh 3 wherein said one or more solid phase supports are 
microparticles and wherein said step of sorting said cDNA molecules onto the 
microparticles produces a subpopulation of loaded microparticles and a subpopulation 

35 of unloaded microparticles. 

5. The method of claim 4 further including a step of separating said loaded 
microparticles from said unloaded microparticles. 



-55- 



wo. 97/13877 



lI^/yS96/16342 



6. The method of claim 5 further including a step of repeating said steps of 
sampling, sorting, and separating until a number of said loaded microparticles is 
accumulated is at least 10,000. 

5 

7. The method of claim 6 wherein said number of loaded microparticles is at 
least 100,000. 

8. The method of claim 7 wherein said number of loaded microparticles is at 
10 least 500,000. i , - ^ ' ' - -) 



:o 



:> 



: } 



^ ' method of claim 5 furtfierjncludihg a:^ep 6f repeating said steps of 

sampling, sorting, and separating until a number of said loaded microparticles is . 
acciunulated is sufficiernt toiestimate^the relatiyerabUndance of a cDfsIA molecule . 
present in said population at a frequency within the range of from 0.1% to 5% with a 
95% confidence limit no larger than 0.1% of said population. . 

1 0. The method of claim 4 wherein said test organism is a mammalian tissue 



- culture. 



1 1 The method of claim 1 0 wherein said mammalian tissue culture comprises 
hepatocytes. ; , 

( 

12. The method of claim 4 wherein said test organism is an animal selected from 
the group consisting of rats, mice, hamsters, guinea pigs, rabbits, cats, dogs, pigs, and 
monkeys. 



13. The method of claim 12 wherein said one or more tissues are selected from the 
group consisting of liver, kidney, brain, cardiovascular, thyroid, spleen, adrenal, large 
0 intestine, small intestine, pancrease urinary' bladder, stomach, ovary, testes, and 
mesenteric lymph nodes. 



14. A method of identifying genes which are differentially expressed in a selected 
tissue of a test animal after treatment with a compound, the method comprising the 
steps of: 

administering the compound to a test animal; 
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. ^ extracting a population of mRNA molecules from the selected tissue pf the. , 



test animal; 

- :forming a population of cDNA molecules from the population of mRNA 
molecules such thai each cDNA molecule has an oligonucleotide tag attached, the 



, it*. 



•i. • 



5 oligonucleotide; tags cross-hybridizing set; ^ ' 

. , . sampling the population of cDNA molecules such that substantially all 

. , ... different cDNA molecules have different oligonucleotide tags attached; _ 

sorting the cDN A molecules by specifically hybridizing the bligohucleotide 
tags with their respective complements, the respective complements being attached as 
0 ' 'uniform pbpulationis of substantially identical complements in spatially discrete^ - ^ ■ 
' 'regions oh one or more solid phase supports; 



determining the nucleotide sequence of a portion of each of the sorted cDNA 
molecules to form' a fiTSqubncy distribution of expressed genes; arid 

identifying genes expressed in response to administering ^e compoimd by 

comparing the frequencing distribution of expressed genes of the selected tissue of the 

' ' . > -It. ' , t » -, ' ' , i f ■ \ • 

test animal with a frequency distribution of expressed genes of the .selected tissue of a 
control animal. ^ K 



- vl5. The method of claim 14 wherein said oligonucleotide tag and said 
0 complement of said oligonucleotide tag are single stranded. ' 

.16. The method of claiin 15 WHerein said oligdnucleotide^u^ 
plurality of subunits, each subunit consisting of an oligonucleotide of 3 to 9 
nucleotides in length and each subunit being selected from the same minimally cross- 
5 hybridizing set. 

1 7. The method of claim 16 wherein said one or more solid phase supports are 
microparticles and wherein said step of sorting said cDNA molecules onto the 
microparticles produces a subpopulation of loaded microparticles and a subpopulation 

0 of unloaded microparticles. 

18. The method of claim 17 frirther including a step of separating said loaded. 

microparticles'from said unloaded microparticles. 

- * * . i * 

• * ■ • \ ' ■ ■ 

- - * _ '*•,.■. . !» i^' 

5 19. The method of claim 1 8 further including a step of repeating said steps of 
. sampling, sorting, and separating imtil a number of said loaded microparticles is 
accumulated is at least 10,000. 



-57- 



22; ^ J, The method of claim 1 8 further including a step of repeating said steps of 
sampling, sorting, and separating until a number of said loaded microparticles is 
accumulated is sufficient to estimate the relative abundance of a cDNA molecule 
present in said population at a frequency within the range of from 0. 1% to! 5% with a 
95% confidence limit no larger than 0. 1 % of said population. 

23. The method of claim 1 7 wherein said test animal is selected from the group 
consisting of rats, mice, hamsters, guinea pigs, rabbits, cats, dogs, pigs, and monkeys. 

24. The method of claim 23 wherein said selected tissue is selected from the 
group consisting of liver, kidney, brain, cardiovascular, thyroid, spleen, adrenal, large 
intestine, small intestine, pancrease urinary bladder, stomach, ovary, testes, and 

- mesenteric lymph nodes. 

25. A use of the technique of massively parallel signature sequencing to determine 
the toxicity of a compound in a test organism, the use comprising the steps of: 

administering the compound to a test organism; 

extracting a population of mRNA molecules from each of one or more tissues 
of the test organism and forming a population of cDNA molecules for each of the one 
or more tissues; 

determining the nucleotide sequence of a portion of each of the cDNA 
molecules of each separate population using massively parallel signature sequencing 
to form a frequency distribution of expressed genes for each of the one or more 
tissues; and 

correlating the frequency distribution of expressed genes in each of the one or 
more tissues with the toxicity of the compound. 

26. The use of claim 25 wherein said test organism is a mammalian tissue culture. 



27. The use of claim 26 wherein said mammalian tissue culture comprises 
hepatocytes. 
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28. The use of claim 25 wherein said test organism is an animal selected from the 
group consisting of rats, mice, hamsters, guinea pigs, rabbits, cats, dogs, pigs, and 
monkeys. 

5 29. The use of claim 28 wherein said one or more tissues are selected from the 
group consisting of liver, kidney, brain, cardiovascular, thyroid, spleen, adrenal, large 
intestine, small intestine, pancrease urinary bladder, stomach, ovary, testes, and 
mesenteric lymph nodes. 

10 30. A use of the technique of massively parallel signature sequencing to identify 
genes which are differentially expressed in a test organism after treatment with a 
compound and which are correlated with toxicity of the compound, the use 
comprising the steps of: 

administering the compound to the test organism; 
1 5 extracting a population of mRNA molecules from a selected tissue of the test 

organism and forming a population of cDNA molecules; 

determining the nucleotide sequence of a portion of each of the cDNA 
molecules using massively parallel signature sequencing to form a frequency 
- distribution of expressed genes; 
20 identifying genes expressed in response to administering the compound by 

comparing the frequencing distribution of expressed genes of the selected tissue of the 
test organism with a frequency distribution of expressed genes of the selected tissue 
of a control organism; and 

determining whether the genes expressed in response to administering the 
25 compound are correlated with toxicity of the compound in the test organism. 



59 



wo 97/13877 



PCT/US96/16342 



1>2 



100 



N 



•J ' 



'"' Generate Table-Mn of 
' all possible subunits of >^ 
desired'length and composition 



Q 



J 1.. . 



110 



i 



Select 'inltiarsubunrt 

'Si(i=i)^'' ' 



120 



Compare subunit Si to 
successive subunits in 
Table Mn from S+1 to 
end of Table 



1 



125 



Saive subunit in 
Table Mn+1 



Does 
subunit meet 
criteria?. 



140 



Is this the 
last subunit 
in Mn? 



Yes 



Discard subunit 
Go to next subunit 



150 



140 



Replace Mn 
with Mn+1 




Is this the 
last subunit 
in Mn?. 



No 



1 



-.if. • - 1 ■ 



' 1 



SUBSTITUTE SHEET (RULE 26) 



wo 97/13877 



PCT/US96/16342 



2/2 



318 



350 




352 




354 




356 










330 




332 




334 




336 


>^ 


>^ 


>^ 


1 





COMPUTER 



/ 
304 

312 
J— 



UGHT 
SOURCE 




Fig. 2 



SUBSTITUTE SHEET (RULE 26) 



INTERNATIONAL SEARCH REPORT 



IntenMtional tppUeadon No. 
PCT/US96/16342 



A, CLASSinCATION OF SUBJECT MATTER^ . , , ... >. 
IPC(6) : C12Q 1/68; O07H 21/04 ~ " " . - . = ^ - . _ . . 

US CL : 435/6; 536/24,3 i:j/oJX- (^i'lj' L .;'i>.,i' v ;'^ 

According to Intenutioiul Patent CUuification (IPC) or to both liitioiul cUaiiBcation and IPC 



B. nELDS SEARCHED 



Minunum documentition learcbed (clauification tyitem Tblkyw^ by clutiOcation symbob) 



U.S. : 435/6; 536/24.3 



Documentation icaicbed other than minimum documentation to the extent that such documentj are included in the fields searched 



I. . 



' 1 ■ -v 



Electronic dau base consulted during the intemationa] search (name of dau bue andv where practicable, search terms used) 
APS, MEDUNE, BIOSIS, CAPLUS, SOSEARCH 

search terms: Martin, David W., toxic?, differential?, express?, cDNA, nnRNA, RNA. oene#, hybrid?. 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category* 



Citation of document; with ^indication, where a p prop ria te; of the relevant 'passages 



Relevant to claim No. 



CHETVERIN et al. Oligonucleotide arrays: New concepts and 
possibilities. Bio/Technology. 12 November 1994; Vol. "12/ 
pages 1093-1099, especially pages 1095-1096. 

BRENNER et al. Encoded combinatorial !,cheniistry. 
Proceedings of the National Academy of Sciences USA. June 
1992, Vol. 89, pages 5381-5383. 



i-30 



1-30 



MATSUBARA et a|. cDNA analyses in the human genome 
project. Gene. 15 December 1993, Vol. 135, No. 1-2, pages 
265-274. 



1-30 



' 1 



I x| Fulther documenU are listed in the continuation of Box C. | | See patent family 



annex. 



Speckl c t i cf orica of 



'V 



of 4ie «t which ii not oooulcred 

lo be of paftienkr 

cvUer donwiml publiihcd oo or afker the mit niiliriMl fUmf dsle 

drwimrm which owy iuow 6ovtm oa phority chun(i) or which u 
ciled to ftiMiih dw pnhbeatioo dale of isoidker dmtioo or olher 
(MipectM) 



icfcrrvg to n onl diadoaure. use. exhihiboo or 

filiBf date but hter than 



docmcBt pubUahcd prior lo ttie 
Ibe pr>oniy date 



"T later dorummi publi«faal after the mtenwtional fUiag date or prtority 

date aod ikA in coafUct with the applicaikm but cited lo VDdefstaod the 
priacipk or theory underlyBi the BVCBlioa 

*X* docuBGBl of particular relevaace; the claimed invcatioa canool be 

cooaiderad oovel or canoot be cooaiderad to iovolve an inventive atep 
when d» docuneot ia taken alooc 

'Y* document of paiticular reievaoce: the churned invcntioD canool foe 

cooaiderad lo irvolve an nvenxive atep when the documeal ti 
cofnbbed with ooe or more other auch docuneota. such combinBtioa 
1 obvioua to a pcraoo akiUed io the art 

r of the aame patent family 



Date of the- actual completion of the international search 



27 JANUARY 1997 



Name and mailing addreat of the ISAAJS 
Conimistioner of Piioou and Trademarfca 

Box PCT 

Washtngtoo, D.C. 20231 
Facsimye So. (703) 305-3230 



Date of mailing of the international search report 

19 FEB 1997 



IT 



Authorized offtce. 

1 

SCOTT D. PRIESE 



* r 



( 



Telephone No. (703) 308-0196 



Form PCT/ISA/2tO (second shee()(July 1992)* 



INTERNATIONAL SEARCH REPORT 



C (ConUnuatm).,. DOCUMENTS CONSIDERED TO BE RELEVANT 



Imemitioiul application No. 
PCTAJS96/16342 



Ctfepwy * jK'iC itrtiba of itociMncal.:wiih ndkatktn,: where ^ipnpmle,. of the relevuu puugei 

;W0 95/21944, Al (SMTTHRLINEiBEECHAM. CORPORATION) , < 1 1-30 
17 August 1995, page 4, lincs.l-4,,page,5, lines 31-37, jage 17, , 
lines 15-27, page 18, lines 30-i35, page 20, line 23 to page 21, 

'line ''''''' " i-- • ' ^' ' 



' i^clevim to claim No. 



.J. 



•.lit i ^. » 



.,■3' 



-1 



< - 1 



! J. 



(1 r 



Form PCT/ISA/210 (continuation of second thect)(July 1992}* 



